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I. INTRODUCTION 



A. BACKGROUND 

As costs for material and energy continue to rise, the 
need for smaller, more efficient heat exchangers continues 
to grow. For those applications involving marine vehicles, 
size and weight limitations dictate the use of small, 
highly-efficient heat exchangers. Reducing the size and 
weight of steam condensers used aboard 0. S. Navy ships for 
propulsion and electrical power generation would result in 
lower material costs, and help to aleviate the cramped 
conditions so typical of machinery spaces. Thus, there is a 
strong motivation for continued study in this area. 

The ef fectiveness of condensers is limited by the water- 
side, vapor-side, and wall thermal resistances of the 
condenser tubes. Reducing any one of these will contribute 
to increased hea t- transfer performance, and smaller physical 
size of condensers. Methods of enhancing the vapor-side 
coefficient include the use of "roped" tubes, fluted tubes, 
drainage strips attached on the tubes, finned tubes, and 
coatings applied to enhance dropwise condensation. This 
thesis concentrates soley on finned tubes. 

Since the late 1940s, externally-finned tubes have been 
used to increase the vapor-side heat-transfer coefficient of 
tubes used in refrigeration systems; but condensers used in 
steam systems, such as shipboard propulsion plants, continue 
to use smooth tubes. The high surface tension of water, 
which leads to its tendency to flood the area between fins, 
has resulted in a widely-held belief that such tubes are 
inappropriate for use in steam systems. Recent studies 
[1,2], however, have shown that finned tubes can signifi- 
cantly enhance the heat-transfer rates in such systems. 
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The theoretical treatment of condensation on finned 
tubes is extremely complex due to the large number of vari- 
ables and physical mechanisms involved. The interaction of 
gravitational and surface-tension forces lead to complex 
three-dimensional flow patterns, which are further dependent 
on fin spacing, height, and thickness. Other variables 
include heat flux, vapor shear, tube diameter, fin shape and 
fluid properties just to name a few. In view of the above, 
any theoretical models will require numerous simplifying 
assumptions, and require complex computer solutions 
involving implicit numerical techniques [3]* To confirm the 
validity of theoretical models, reliable experimental data 
which cover a wide range of relevant parameters must be 
obtained. The availability of such data may then lead to 
the development of fairly simple experimental correlations, 
which could be used in the design stage to predict the heat- 
transfer performance of finned tubes. 

This thesis effort is a continuation of research being 
conducted at the Naval Postgraduate School (UPS) under a 
grant from the National Science Foundation. The basic test 
apparatus used to collect experimental data was built by 
Krohn [4], Graber [5] provided the majority of instrumenta- 
tion, and took preliminary data as the system experienced 
problems with non-condensing gases and partial dropvise 
condensation on copper tubes. Poole [6] made further 
improvements on the apparatus as well as on the instrumenta- 
tion and, most importantly, assured a leak-free apparatus. 
Unfortunately, he did not have sufficient time to produce 
useful data, mostly due to the considerable time spent in 
systematically locating and fixing leaks, and due to the 
partial dropwise condensation problem that had not been 
solved. Using this system, Georgiadis [7] was finally able 
to obtain complete filmwise condensation on copper tubes. 
He took data on a number of finned tubes with fins of 
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rectangular shape, as well as on smooth tubes. The repeat- 
ability of data obtained by Georgiadis demonstrated the 
accuracy of the test apparatus and associated instrumenta- 
tion which was used essentially without modification for 
this investigation. 

The overall objectives of the present program at NPS 
includes the testing of: (a) tubes with rectangularly- 

shaped fins to find the best fin spacing, thickness, and 
height, (b) tubes with various fin shapes to find the best 
geometry that will maximize heat transfer, (c) tubes with 
different tube-metal thermal conductivity, and (d) the 
effect of vapor shear on finned tubes. Georgiadis [7] 
tested a total of 25 tubes, which were divided into five fin 
spacings (0.5, 1.0, 1.5, 2.0 and 4.0 mm), five fin thick- 
nesses (0.5, 0.75, 1.0, 1.25, and 1.5 mm), and two fin 

heights (1.0 and 2.0 mm). At the conclusion of Georgiadis' 
test program, eight new tubes with fin heights of 0. 5 and 
1.5 mm, fin spacings of 1.0, 1.5, 2.0, and 4.0 mm, and a 

fin thickness of 1.0 mm, remained to be tested to complete 
the sequence of tubes with rectangularly-shaped fins. The 
testing of these tubes and additional tubes toward objec- 
tives (b) , (c) , and (d) were the primary goals of this 

thesis effort as listed in the next section. 

B. OBJECTIVES 

The main objectives of this thesis are as follow: 

1. Take data on tubes with rectangularly-shaped fins of 
various fin spacings and fin heights to augment 
previous data [7 ], 

2. Take data on tubes with fins of different shapes 
(triangular, trapezoidal, parabolic, etc.), 

3. Take data on commercially-available finned tubes. 
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4. Take data on the "optimum," rectangularly-shaped 
finned tube at different vapor velocities, 

5. Take data on a stainless steel tube with "optimum," 
rectangularly -shaped fins, and 

6. Develop a preliminary correlation based on data for 
tubes with rectangularly-shaped fins. 
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II. THEORETICAL TREATMENT OF CONDENSATION ON HORIZONTAL 

FINNED TUBES 



A. FILM CONDENSATION 

When filmwise condensation of a vapor tak. 
smooth tubes, a thin layer of condensate is 
thickens with increasing distance around the 
the tube. This condensate layer creates a th 
tance, which can limit the heat- tra nsf er perfor 
tube. This film thickness and its thermal resis 
reduced by the use of external, radial fins whir 
tion to a surface area increase, promote sur 
effects. 

In 1984, Yau et al. [1] measured the 
provided by copper finned tubes over smooth tub 
wise condensation of steam. Similar exp 
Wanniarachchi et al. [2] also in 1984 confirm 
observed enhancements were greater than could 
by the increased surface area alone. This 

enhancement may be a result of the surface-te 
which act to thin the condensate film. Figure 2 
cally depicts this phenomenon. 

The effect of surface tension on pressure a 
face between a liguid and vapor is inversely pro 
the radius of curvature of the interface. If a 
two radii of curvature at right angles (i.e., r 
can le shown that: 



AP = o 
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formed which 
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ermal resis- 
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viier e 



o - surface tension of condensate, and 
AT = the pressure difference. 



fin 



Condensate 




Figure 2.1 Schematic of Condensate Profile at Upper 

Tube Surface. 

For the case of a finned tube, the radius of curvature 
around the outer fin perimeter is very large compared to the 
radius of curvature around the fin cross-section profile 
(Figure 2.1). Feferriny bach to equation (2.1), it can be 
seen that the smaller radius of curvature term will domi- 
nate, so surface-tension effects around the fin perimeter 
may be neglected. 

Because of the convex shape of the condensate film at 
point A, the pressure within the film at this point is 
greater than the surrounding vapor pressure. In a similar 
manner, the pressure at point B is less than the surrounding 
vapor pressure owing to the concave shape of the condensate 
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film. The relatively flat shape of the condensate film at 
point C leads to a pressure essentially equal to the 
surrounding vapor pressure. These pressures are given by: 
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(2.4) 



where 



P v = surrounding vapor pressure, 

P A , P s' p c = liquid pressure at points A, B, and C, and 

r A , r g, r c = radius of curvature of the condensate 
film at points A, B, and C. 



At point B, the radius of curvature is small, so the pres- 
sure at point B is less than the pressure at point C. 
Further, from equations (2.2) and (2.3), the pressure at 
point A is greater than at point B. In reality, the pres- 
sure gradient within the condensate film varies along the 
height of the fin due to the continuously varying radius of 
curvature from the fin tip down to the base [3]; however, to 
simplify the treatment of condensate flow toward the fin 
base, the pressure differences between points A and B and 
points C and B may be simply written as: 



22 



AP 



AT 



o 




( 2 . 5 ) 




o 



( 2 . 6 ) 



where 



Ap ab , AP cb = prpssure difference between points A and B, 



As can be seen, these pressure differences are positive. 



relatively large mass of condensate at point B, gravita- 
tional forces dominate over surface-tension forces causing a 
flow of condensate around the tube perimeter at the fin 
base. In this manner, a condensate run-off channel is 
formed at point B, and the improved drainage of condensate 
thins the film between fins and on the fin surface. This 
thinning, in turn, reduces the thermal resistance through 
the condensate film, thus producing an enhancement in addi- 
tion to the gain in surface area. 

It should be stressed at this point that the shape of 
the condensate film is highly dependent on the fin geometry. 
Clearly, then, the fin shape must be taken into considera- 
tion in order to maximize the beneficial effects of surface 
tension. 

The above gains may be partially or totally offset, 
however, by the tendency of condensate to flood the area 
between fins, especially on the lower part of the tubes. 



and points C and B. 



resulting in condensate flow toward point B. Due to the 
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The extent to which the film floods the tube is defined by 
the condensate retention angle ( ij> ) , as shown in Figure 2.2- 
In the flooded portion of the tube, the relatively thick 
condensate film increases the thermal resistance, thereby 
leading to a degradation in the heat-transf er performance. 
The condensate retention angle is highly dependent on fin 
spacing, so very closely-spaced fins may lead to completely 
flooded tubes. 

A number of studies have been conducted to investigate 
the performance of finned tubes, and these are presented in 
the next section. 

B. CONDENSATE RETENTION 

In 1946, the first measurements of condensate retention 
were made by Katz et al. [8]. These measurements were made 
under static conditions (i.e., no condensation taking 
place), using water, aniline, acetone, and carbon tetrachlo- 
ride on a number of tubes with different fin densities (276 
to 984 fins/m), and fin heights (1.2 to 5.7 mm). It was 
shown that as much as 1005?. of the tube surface could be 
flooded by retained condensate, depending mainly on the 
ratio of surface tension to liquid density and on the fin 
spacing. 

Eight years later, Gregorig [9] recognized that surface- 
tension effects could play an important role on a vertical, 
fluted surface which uses a minimum radius of curvature at 
the flute tip, which gradually increases toward the trough. 
As discussed in section A above, this variable radius of 
curvature results in a condensate film pressure that 
decreases toward the trough. In this manner, the condensate 
layer on the convex region of the flute is thinned signifi- 
cantly, thus pushing the condensate into the troughs (Figure 
2.3) where gravitational forces result in enhanced drainage. 
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External Diameter of fins 



Root Diameter of fins 



7 ube Wa 




Retained Condensate 



Rentention Angle (C) 



Figure 2.2 Schematic of Condensate Retention 
on Finned Tubes. 
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Since the heat- transfer coefficient is inversely propor- 
tional to the thickness of the condensate layer, heat- 
transfer performance at the trough is degraded, but to a 
lesser degree than the enhancement experienced by the crest. 
The resulting average heat- transfer coefficient, therefore, 
shows a significant enhancement. 



THIN FILM REGION 




Figure 2.3 Cross Section of the Gregorig Surface 

In 1381, Rudy and Webb [10] made measurements of conden- 
sate retention angles on finned tubes with three different 
fin densities (748, 1024, 1378 f ins/ieter) using water, 
E- 1 1 , ar.d n-pentane as the working fluid. later, in 1583, 
Pudy and Webb [11] developed an analytical model to predict 
the fraction of tube surface that is flooded during conden- 
sation on a horizontal, integral-fin tube. They found that 
the vertical-rise height of condensate on a finned tube was 
the same as that obtained on a vertically oriented, flat 
finned plate that was obtained by splitting and unrolling an 
identical finned tube. Based on this observation, they used 
capillary equations that predict the liquid rise on a 
vertical U-shaped channel, and assuming negligible vapor 
shear, they predicted the condensate retention angle to be: 
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i|' = cor. 



(2.7) 



-1 



? o (2 e — t) 
Pf r “ A n o ~ 



wher e 



e = fin heig lit , 
t = fin thickness, 
p f = density of condensate, 
g = acceleration of gravity, 
s = fin spacing, and 
D 0 = root diameter of tube. 



This result shows that the condensate retention angle 
increases both with increasing fin density, and with 
increasing surface tension-to-density ratio. It should be 
kept in mind, however, that the model was based on a 

vertical surface, so it should only be used for angles less 
than about 30 degrees. For angles below this, experimental 
results involving the use of water, R-11, P-12, ammonia, and 
n-pentane were predicted to within 10 percent. It is worth 
mentioning that equation (2.7) had been derived in 1982, 
according to Russian literature, by Rifert [ 12]. 

Like Rudy and Webb [10], Owen et al. [13] also recog- 
nized the need to consider condensate retention while 
analyzing finned tubes. In order to correct for this, an 
assumption was made that the condensate retention angle was 
independent of condensation rate, so a static analysis was 
performed. A simple force balance between surface tension 
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arid gravitational forces resulted in an equation for the 
condensate retention angle as shown below: 



ij' = ros 



-1 



4 a 

Pf n « r>f 



( 2 . 8 ) 



where 



D f = overall diameter of fins. 



This equation is the same as equation (2.7), except that 
equation (2.8) is independent of fin thickness ( t ) . 

In 1983, Honda et al. [3] performed a theoretical anal- 
ysis to determine the condensate retention angle. Using an 
iterative numerical scheme, they found the solution to agree 
with equation (2.8) obtained earlier by Owen et al.. Using 
data of their own, Honda et al. verified a close agreement 
between the predicted and experimental values of the conden- 
sate retention angle. 

Yau et al. [ 1 ] conducted experiments to determine the 
effects of removing retained condensate from finned tubes by 
installing thin metal drainage strips attached edgewise to 
the bottom tube surface. Experiments were conducted where 
the vapor velocity and fin pitch were varied on finned tubes 
with and without drainage strips. They showed a significant 
reduction in the in the condensate retention angle when the 
finned tubes were fitted with drainage strips. Condensate 
retention angles for condensation of steam, ethylene glycol, 
and Ii-113 on finned tubes with drainage strips were found to 
fit the empirical relation listed below to within ± 10 
percent : 
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^ = cos 



In 1985, Rudy and Webb [14] expanded their 1983 model in 
order to predict the condensate retention angle for hori- 
zontal finned tubes with fins of arbitrary shape. As 
before, their model was based on capillary equations that 
predict the amount of liquid rise, and negligible vapor 
shear. The resulting equation for the condensate retention 
angle is given by: 



1 - 



1 .66 o 
Pf c r, n f 



(2.9) 



if* = cos 



l - 2 o (P £ - to 

Pf 4 Kt h + si e - Apl 



( 2 . 10 ) 



where 



P^ = wetted perimeter of fin cross section, 
t^ = fin base thickness, and 

Ap = profile area of fin over fin cross section. 



In order to test the validity of equation (2.10), it was 
compared to experimental data for horizontal finned tubes 
with fins of primarily trapezoidal cross section, and fin 
densities ranging from 630 to 1614 fins/m. Steam, K — 1 1 , and 
n-pentane were used as the condensing fluids. Equation 
(2.10) was shown to predict the experimental data to within 
± ID percent. Note that equation (2.10) reduces to equation 
(2.8) for fins of rectangular cross section. 



29 



C. THEORETICAL MODELS 



In 1947, Beatty and Katz [15] performed a number of 
experiments where various refrigerants were condensed on 
single finned tubes. To predict the film coefficients, they 
started with the Nusselt equations for condensation on a 
horizontal tube and on a vertical plate, and with the 
respective proportions of horizontal tube area (between 
fins) and vertical fin area, they calculated the total heat 
transfer rate. From this, an average heat-transfer coeffi- 
cient was obtained based on an equivalent tube diameter. 
Their final expression is simply the Nusselt equation for a 
smooth horizontal tube, but with the tube diameter replaced 
with the equivalent tube diameter. The leading coefficient 
was modified to fit their experimental data as shown below: 



h = n.6«n 
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A cff = A s + n A f 



(2. 13) 



2 2 

ti rn f - n n ) 

4“D 7 (2. 14) 



where 



A Q = surface area of smooth tube, 
h = average vapor-side heat-transfer coefficient. 
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k 



f 



p v 



h 



fg 



v 



f 



AT 



D 



e 



n 



A 



s 



A 



f 



A eff 



thermal conductivity of condensate, 
density of vapor, 

specific enthalpy of vaporization, 

viscosity of condensate, 

vapor-side temperature drop, 

equivalent tube diameter, 

fin efficiency, 

surface area of smooth tube, 

total surface area of finned tube, and 

effective area of finned tube. 



Accuracy of this equation was claimed to be better than 11 
percent for a wide variety of nonagueous fluids. It should 
be noted, however, that this model ignores the effect of 
surface-tension forces; so, it is valid only for low- 
surface-tension fluids such as refrigerants. 

Analytical and experimental studies of condensation on 
horizontal tubes with tra pezoidally-shaped fins were 
performed by Zozulya, Karkhu, and Borovkov [16,17] some 
years later in the 1970s. Their experiments confirmed the 
need to consider surface-tension forces in addition to 
gravitational forces when developing mathematical models. 
The analytical solutions were based on the following assump- 
tions; 1) the thin condensate film on the fins was treated 
as a laminar boundary layer with a pressure gradient along 
the fin profile caused by surface-tension forces; 2) the 
effect of gravitational and inertial forces on the motion of 
the film along the side surfaces of the fins into the 
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condensate-filled trough was neglected; 3) the motion of 
condensate in the trough area is laminar and produced by 
gravity; and 4) no condensation takes place on the flooded 
portion of the tube. Differential equations were obtained 
for the height of condensate between fin bases, and tempera- 
ture distribution along the fin height. With the aid of a 
computer, solutions for these equations were obtained using 
numerical methods, leading to the expressions below; 



h 



F r AT 



(2. 15) 
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JL + b + 
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(2. 16) 



G 



2 4 3 

pf g e sin (tt-i}/) 

TTTf 



7 + 



2 G 



tan( a) 



(2. 17) 



where 

G = condensate flow rate, 

t t = fin tip width, 

a = fin semivertex angle, and 

7 = dimensionless depth of condensate 
between fins. 



Although the above equations appear to be fairly simple, it 
must be kept in mind that rather complex computer solutions 
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are required to obtain values for Z. To check the validity 
of equation (2.15), (2.16) and (2.17), experiments were 
performed with four different finned tubes to condense both 
steam and R— 113 vapor. Agreement to within 5 percent was 
reported. 

In 1S77, Nader [18] presented an analytical model to 
predict condensation on a vertical, single fin attached to a 
horizontal tube. This model differed from those discussed 
above in that the fin temperature was allowed to vary along 
the height of the fin. For the special case where the fin 
temperature is constant, the Nusselt eguation was obtained. 
A "condensation efficiency" was then proposed to account for 
a variable fin temperature, and was defined as: 



<5 = 






(2. 18) 
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1 kf Pf AT 



(2. 19) 



F = C km 
? ? e k f 



( 2 . 20 ) 



where 



C 1 = constant of proportionality, and 

k = thermal conductivity of fin. 
m 



Pith this definition, the rate of condensate formation and 
the rate of heat transfer (with variable fin temperature) 
were found to be: 
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1/4 



( 2 . 21 ) 
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Q = 1.8856 k^ F <I>^ AT 



( 2 . 22 ) 



where 



G f = rate of condensate formation, and 

Q = total heat-t ransfer rate. 

Thus, with a table of values for and F 2 obtained numeri- 
cally, the rate of heat transfer and condensate formation 
could be calculated by a fairly simple expression provided 
the condensation efficiency is greater than 0.8. However, 
for values below 0.8, the accuracy of equations (2.21) and 
(2.22) diminishes quite rapidly. 

In 1979, Patankar and Sparrow [19] treated the case of a 
fin attached to a cooled vertical plate or a cylinder as a 
three-dimensional problem. A thin fin was assumed, and 
temperature variations across the thickness were ignored as 
was lengthwise conduction. Additionally, all heat transfer 
was assumed due to condensation only. Compared to the ideal 
fin solution, results obtained by Patankar and Sparrow indi- 
cated that a significant error existed when using the ideal 
fin solution on all but very long fins. Based on their 
results, the following relations were proposed: 
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where 



g = heat flux, and 
z = axial coordinate. 



According to Patankar and Sparrow, these equations should be 
valid for most practical applications; however, no experi- 
mental data were available at the time to confirm this. 

In 1980, Rifert [20] analyzed condensation of vapor on 
horizontal finned tubes enhanced by the effect of surface- 
tension forces. In his analysis, he divided the tube into 
flooded and unflooded zones, and solved a two-dimensional 
form of the energy equation for each zone. Ihe mean heat 
flux was then determined by integrating over each zone and 
the tube length. In cases where condensate is retained in 
more than half of the tube perimeter, Rifert points out that 
a three-dimensional form of the energy equation must be 
used. Solutions to these equations revealed that, in most 
cases, the fin temperature is very nonuniform, so the mean 
integral heat-transf er coefficient and its correlations 
should not be used. The above-mentioned zone-by-zone anal- 
ysis was, therefore, recommended. 

Based on their study of condensate retention mentioned 
earlier in section B, Rudy and Webb [10] proposed that the 
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Beatty and Katz [15] model be modified to account for 
condensate retention, as shown below: 



h 



= b 



PK 




(2.25) 



where h is computed using equation (2.11). This equation 
neglects any heat transfer through the flooded portion of 
the tube. As a result, it was shown to underpredict the 
average heat-transfer coefficient of condensing R-11 by as 
much as 30 percent when a significant amount of retained 
condensate was present. 

Using the condensate retention angle they had developed 
earlier (see section B) , Owen et al. [13] divided a hori- 
zontal finned tube into an upper unflooded section, and a- 
lower flooded section. Unlike Rudy and Uebb [10], however, 
they accounted for heat-transfer through the flooded portion 
of the tube by noting that condensation occurs on both the 
surface of the retained condensate and the fin tips. An 
effective thermal resistance was obtained over this portion 
of the tube by assuming parallel paths for heat- transfer 
through the fins and retained condensate. The respective 
thermal resistances, which were added in parallel, were then 
added in series with the tube wall thermal resistance. The 
resulting heat-transfer coefficient was then found to be: 



h 




TT 



h 

U 




(2. 26) 
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where h u is computed using equation (2.11), and 




k f 3 p f (p f - p v i r h f p 




(2. 27) 



k _ r = effective thermal conductivity 
e 1 1 “ 

This equation was found to precict all the available data 
with an accuracy of better than 30 percent. 

Rudy and tlehb [21] used equation (2.7) a short time 
later while developing an equation for the heat- transfer 
coefficient. The model for this undertaking was based on 
sur f a ce-tension- dr iven radial flow of condensate on the fin 
surface, with a linear pressure gradient along the fin 
height, and gravity-drained flow of condensate in the 
channel between fins. The Nusselt equation for horizontal 
tubes was used for the tube area between fins, while the fin 
surface was treated by an equivalent gravity model developed 
by Rebb et al. [22] and Rudy [23] earlier. As before, vapor 
shear and heat transfer through the flooded portion were 
neglected, yielding a heat-transfer coefficient of: 
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where 



I = length of tube, 

A bt = surface area of tube between fins, and 
A ft = fin surface area. 

This expression provided an accuracy of better than 1055 for 
condensation of R-11 on short, finely-spaced fins, but the 
accuracy dropped sharply for larger fins spaced further 
apart. This was apparently due to the assumed linear pres- 
sure gradient on the fin surface, and breakdown of the model 
as gravity forces became dominant. The use of equation 
(2.28) should be restricted to fin densities from 1200 to 
1400 fins/m, and fin heights of less than 1 mm. 

In 1984, Honda et al. [24] developed a model to predict 
the heat-transfer coefficient that took surface-tension 
effects and the non- isothermal behavior of fins into 

account. In order to obtain a solution, however, iterative 
finite difference techniques that required a considerable 
amount of computer time were used. Nonetheless, they showed 
that this model was able to predict the average heat- 
transfer coefficient for most available experimental data 
(which included 11 fluids, and 22 tubes) within 20 percent. 
Predictions involving the condensation of steam provided the 
largest errors (as much as 4055) , due in part to the high 
surface tension of water. This reflects the fact that 
surface-tension effects are still not completely accounted 
for. Compared to the other models discussed thus far, 
however, the model of Honda et al. shows considerable 
promise in spite of its complexity. 

In 1983, Adamek [27] defined a family of condensate 
surface profiles whose curvatures are given by: 
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The parameter £ in eguation (2.29) char act erizes the aspect 
ratio of the fin cross section. A number of £ values and 
their corresponding condensate surface profiles are shown in 
Figure 2.3. With the above information, Adamek obtained an 
expression for the heat-transfer coefficient of the fin 
surface as shown below: 

1/4 

0 h fp, °m S n Pf (5 + 1) ] /0 

Uf k f AT TT+ 2) 7 U.JU) 

These fin profiles were used by Webb et al. [25] and Rudy 
et al. [26] in 1985 to develop a model for condensation on 
horizontal finned tubes. This model divided the tube into 
flooded, and unflooded regions, and further divided the 
unflooded region into tube (area between fins) and fin 
areas. To predict heat transfer from the unflooded fin 
surfaces, they used the heat-transf er coefficient proposed 
by Adamek. The Adamek fin profile used to approximate their 
trapezoidal fins is shown in Figure 2.4. The Nusselt 



h r = 2.140 



'm 
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Figure 2.3 Adamek [27] Condensate Surface Profiles 




Figure 2.4 Fin Geometry for the Webb et al. Hodel [22] 
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equation, written in terms of the film Reynolds number, was 
used to predict heat transfer from the unflooded tube 
surface between fins. Writing the Nusselt equation this way 
takes into account the additional condensate thickness that 
results from the drainage of condensate formed on the fin 
surface. For the flooded region, a two-dimensional computer 
code was used to solve for the heat flux into the tube-side 
coolant ( q b2 ) , which was compared to the limiting case of 
heat flux (q bl ) assuming zero fin thickness. This ratio 
(<t> - % 2 A J bl ) was use< 3 along with a linear temperature 

profile across the condensate film to establish a heat- 
transfer coefficient for this region. The resulting average 
heat- transfer coefficient is given below: 
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where h f is calculated using equation (2.30), and 



A , = smooth tube area based on the fin diameter, 
sf 

Re f = film Reynolds number. 
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ratio of tube side heat flux with fins to tube 
side heat flux with fins of zero thickness. 



A detailed step-by-step procedure, recommended for obtaining 
a solution of equation (2.31), is provided in Appendix A. 

To allow for heat transfer from the fin tip, and to 
account for decreased heat transfer due to the thick conden- 
sate film in the trough area, Webb et al. replaced ^ with 
s m + t/2 - 6 r , where 6 r is the average condensate film thick- 
ness in the trough area in the unflooded region: 
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Ill- DESCRIPTION of test apparatus 



A. TEST APPARATUS 

The same test apparatus used by Georgiadis [7] was used 
for this investigation. A schematic of this apparatus is 
shown in Figure 3.1. Steam was generated in the boiler, 
which consisted of a 304.8 mm (12 in.) Pyrex glass section 
containing ten 4000-Watt, 480-Volt Watlow immersion heaters. 
After passing through a 304.8 mm (12 in.) to 152.4 mm 
(6 in.) reducing section, the steam flowed upward through a 
2.44 m (8 ft.) long section of Pyrex glass piping. At this 
point, a 180-degree bend in the piping re-directed the steam 
downward, where after 1.52 m (5 ft.) of Pyrex glass piping, 
it entered the stainless steel test section illustrated by 
Figure 3.2. The test tube was mounted horizontally in the 
test section behind a view port, which allowed visual obser- 
vation of the condensation process. Steam not condensed by 
the test tube continued downward to the auxiliary condenser, 
where it was condensed by two 9.5 mm (3/8 in.) diameter 
water-cooled copper tubes helically coiled to a height of 
457 mm (18 in.). The condensate was then returned to the 
boiler by gravity. 

Filtered tap water was used to cool the test tube on a 
once-through basis. This water was first collected in a 
large sump with a capacity of about 0.4 cubic meters (Figure 
3.3), then pumped through a flow meter, and the test tube by 
two centrifugal pumps connected in series. A valve on the 
discharge side of the second pump allowed the velocity of 
water flowing through the test tube to be varied from 0 to 
4.4 m/s (14.4 ft/sec). The auxiliary condenser was cooled 
by a continuous supply of tap water, which was throttled to 
control the internal pressure of the test apparatus. 
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Figure 3.1 Schematic of Test Apparatus 
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Figure 3.2 Schematic of Test Section (Insert Removed) 



The system used to remove non-condensing gases is shown 
in Figure 3.3. The vacuum pump continually drew a sample 
mixture of any possible non-condensing gases and moisture 
from the auxiliary condenser. This mixture was then passed 
through a heat exchanger, where the moisture was condensed 
before it entered the sump, using the same filtered tap 
water described earlier. 

B. INSTRUMENTATION 

The electrical power input to the boiler immersion 
heaters was controlled by a panel-mounted potentiometer. A 
converter with an input voltage of 440 VAC generated a 
signal which was fed to the data acquisition system to 
calculate the boiler input power. A more detailed descrip- 
tion of the boiler power supply is provided by Poole [6]. 

The internal system pressure was measured manually using 
a tJ-tube, mercury-in-glass manometer graduated in millime- 
ters. Steam, condensate return, and ambient temperatures 
were measured using calibrated copper-con stantan thermocou- 
ples. These thermocouples had an accuracy within ± 0.1 K. 
The cooling water temperature rise was measured by two 
Hewlett-Packard (HP) 2804A quartz thermometers, along with a 
10- junction, series-connec ted copper-constan tan thermopile 
as a backup. Throughout this investigation, the quartz 
thermometers and the thermopile agreed to within ± 0.03 K. 

C. DATA ACQUISITION SYSTEM 

To monitor the system temperatures and boiler input 
power (using the converter signal) , an HP 9826A computer was 
used to control an HP 3497A Data Acquisition System. Raw 
data were processed immediately and, at the same time, 
stored on diskett for reprocessing at a later time. 
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Figure 3.3 Schematic of Vacuum System and Cooling Water Sump 



D. TDBES TESTED 

For this thesis effort, eighteen copper tubes, and a 
stainless steel tube were manufactured. They are presented 
in Table I. The first eight tubes with rectangularly-shaped 
fins were tested in order to complete the sequence of data 
initiated by Georgia dis [7]. This completed the systematic 
variation of fin spacing, thickness, and height so that the 
optimum dimensions could be determined. Figure 3.4 shows a 
photograph of the four tubes with fin thickness of 1.0 mm, 
fin height of 0.5 mm, and spacings of 1.0, 1.5, 2.0, and 4.0 
mm (tube numbers 26 to 29). Also shown is the insert used 
to enhance the inside heat- transfer coefficient, which will 
be discussed in more detail in the next chapter. Figure 
3.5(a) shows a cross section of a rectangular fin with a fin 
height of 1.5 mm, a fin spacing of 1.5 mm, and a fin thick- 
ness of 1.0 mm. In addition to the tubes with 

rectangularly-shaped fins, six tubes were tested to study 
the effect of fin shape on the heat-transfer performance. 
Tubes number 40 and 42 were made of stainless steel to test 
the effect of fin-metal thermal conductivity on the heat- 
transfer performance. Due to limited machine-shop facili- 
ties available, the "parabolic" fins on two of the tubes 
were not truly parabolic, as shown in Figure 3.5(b). A 
cross section of a commercially-available finned tube 
produced by the Wolverine Division of United Oil Products is 
also shown in Figure 3.5(c) for comparison. The smooth 
tubes were tested in order to determine the water-side heat- 
transfer coefficients, and to serve as a comparison for the 
finned tubes. 

E. VACOUB INTEGRITY 

For any condensation experiment, especially when oper- 
ating under vacuum conditions, a leak-free apparatus is of 



48 



Fin Tube 

Tube Fin Fin Base Fin Tip In ter - Fin Fin Root Inside Tube 

Number Type Thickness Thickness Spacing Height Diameter Diameter Material 
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Stainless Steel 
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Figure 3.4 Photograph of Finned Tubes with a Fin 
Height of 0.5 bio and Insert. 




(a) 





(c ) 



Figure 3.5 
Fins, (b) 



Cross-Sectional Photographs of (al Rectangular 
"Parabolic" Fins, and (c) "Wolverine" Fins. 
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vital importance. The presence of any leaks would allow the 
admission of non-condensing gases, whose deleterious effects 
would invalidate any data being collected. To prevent this 
problem, vacuum integrity tests were performed on the test 
apparatus both before and after the data acquisition phase 
of this thesis effort (a time period of about four months) . 
The first test revealed a leak rate represented by a pres- 
sure rise of 1.0 mmHg over a period of five days, while the 
pressure rise during the final test was 1.0 mmHg over a 
24-hour period. These tests were performed at an absolute 
pressure of about 85 mmHg. This negligible leak rate, 
together with continuous venting of the test apparatus 
during data runs, eliminated the deleterious effects of 
non-condensing gases on the heat-transfer results reported 
in this thesis. 
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17. DATA COLLECTION AND BEDUCTION 



A. SYSTEM OPERATION 

In order to ensure filmwise condensation, all tubes were 
treated with a solution containing equal parts of sodium 
hydroxide and ethyl alcohol. The tubes were placed in a 
steam bath and several layers of the solution were applied 
with a tooth brush every 10 minutes. After about one hour, 
the tube surfaces were completely blackened. The same 
procedure was repeated prior to each data run, but the tube 
was only steamed for ter. to fifteen minutes. Following the 
steam bath, the tubes were rinsed with tap water before 
installing them in the test apparatus. Georgiadis [7] 
discusses this procedure in more detail. 

The test apparatus was brought to operating pressure and 
temperature using the procedures provided by Georgiadis [7]. 
Steady-state conditions were assumed to occur once the test 
section steam temperatures, as well as the cooling water 
inlet and outlet temperatures stabilized. Two sets of data 
were then taken for each of the heat fluxes, which were 
determined by the flow rate of cooling water through the 
test tube. Starting at a flow rate of 80 percent (4.44 m/s 
for 19 mm O.D. tubes, and 3.92 m/s for 13.7 and 15.8 mm G.D. 
tubes), flow rates were reduced to 70, 60, 45, 35, 26, and 
20 percent, then brought back up to 55 and 30 percent. 
These settings were selected since they provided nearly 
equally-spaced heat flux values. It should be noted that 
continuous adjustments of cooling water flow through the 
auxiliary condenser were required to maintain system pres- 
sure as flow rates through the test tube were changed. 



53 



Visual observations of condensation taking place on the 
test tube were made on a regular basis to ensure uniform 
filmwise condensation. A more reliable indication of this, 
however, was the trend of data taken as the cooling water 
flow rates were increased back up to 80 percent. If the 
final value of the steam-side heat- transfer coefficient at 
80 percent was significantly higher than those at the begin- 
ning of the data run, dropwise condensation was assumed to 
have occurred, and the data were disregarded. All data 
presented in this thesis displayed less than 555 disagreement 
in the steam-side heat- transfer coefficient between initial 
and final data sets. 

B. THE DROPWISE CONDENSATION PROBLEM 

The purpose of this thesis was to take strictly filmwise 
condensation data, using primarily copper tubes. Due to the 
poor wettability of copper with water, especially when even 
minute amounts of contamination are present, there is a 
tendency to condense steam in the partial dropwise mode. 
Dropwise condensation is far more effective than filmwise 
condensation; so, its presence can lead to large errors in 
the data. 

By following procedures set forth by Georg.iadis [7], the 
dropwise problem was minimized. During the initial testing 
phase of this thesis, however, a small amount of dropwise 
condensation was encountered. This problem was soon traced 
to the boiler water, which had not been changed for some 
time. By changing the water after every three or four runs, 
the problem was eliminated. 

C. STEAM VELOCITY LIMITATIONS 

A major assumption in the Nusselt theory is that vapor 
velocity does not induce shear forces on the vapor-licuid 
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interface. Since system pressure in the test apparatus was 
controlled by condensing steam in the auxiliary condenser, a 
reasonable steam velocity was always present. Since steam 
velocity (i.e. , steam flow rate) depends on boiler power 
input, decreasing the power results in lower vapor veloci- 
ties. Under vacuum conditions, (about 85 mmHg) a maximum 
velocity of about 8 m/s was possible at full boiler power 
(36 kF) . However, at these high velocities, the entire 
apparatus experienced considerable vibration. As discussed 
by Georgiadis [7], a steam velocity of 2 m/s under vacuum 
provided the most stable operating conditions. This corre- 
sponding boiler power (about 9 kW) was low enough to mini- 
mize system vibrations, yet high enough to avoid rapid 
fluctuations in condenser pressure due to the intermittent 
break-up of vapor bubbles in the boiler. To allow a suffi- 
cient amount of steam to reach the auxiliary condenser when 
operating at atmospheric pressure, the maximum boiler power 
was required. This resulted in a steam velocity of approxi- 
mately 1 m/s. 

D. DATA REDUCTION 

All the programs, property functions, and calibration 
data used for this thesis effort were essentialy the same as 
those used by Georgiadis [7] and Poole [6]. Because six of 
the tubes (tube numbers 39 thru 44) tested during this 
effort had inside diameters different than any tube tested 
previously, and two (tubes number 40, and 42) were manufac- 
tured from stainless steel, some minor modifications were 
necessary for the data reduction program (program DRP6) . 
These modifications included options for different tube 
diameters, tube inserts, thermal conductivity, mixing- 
chamber calibrations, and leading constants for the inside 
heat-transfer coefficients. 



55 



The mixing- chamber (see Figure 3.2) calibration was 
required to account for the temperature rise that occurred 
across the mixing chamber, due primarily to turbulent 
viscous dissipation. The purpose of the mixing chamber is 
to ensure a uniform temperature at the cooling water 
discharge, so that accurate temperature measurements can be 
made. Since tubes with different inside diameters, 
requiring different inserts, result in different pressure 
drops, calibrations had to be performed for each inside 
diameter. 

A modified Wilson Plot program [7] was used to compute 
the leading constants for the Sieder-Tate equation used to 
determine the water-side heat-transfer coefficient. Data 
taken on smooth tubes were used to determine the leading 
constants, which were then used to reduce data taken on 
similar tubes with identical internal dimensions. 
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V. RESULTS and discussion 



A. INTRODUCTION 

During this thesis effort, a number of data runs were 
made using the procedures described in Chapter IV. For both 
vacuum and atmospheric conditions, every tube was tested at 
least three times, on different days, to ensure the data 
were repeatable. Complete filmwise condensation was main- 
tained, and the computed non-condensing gas concentration 
was held below 1.0 percent for all accepted data runs. It 
was estimated that an error of ± 0.2 K in the steam tempera- 
ture, or an error of ± 1.0 mmHg in the system pressure would 
lead to an error of ± 1.2 percent in the non-condensing gas 
concentration. This shows that the computed non-condensing 
gas concentration was zero to within the accuracy of temper- 
ature and pressure measurements. As discussed earlier in 
Chapter III, the test apparatus would allow only a negli- 
gible amount of non-condensing gas to be present in the 
apparatus. The non-condensing gas concentration was 
computed for the system temperature and pressure during 
every data run to ensure no major leaks developed. A 
summary of tubes tested between Georgiadis [ 7 ] and this 
thesis effort and the resulting enhancements are provided in 
Table II. 

B. NATER-SIDE HEAT-TRANSFER COEFFICIENTS 

The water-side heat-transfer coefficients, enhanced by 
the use of an insert, were obtained by Georgiadis [7] for 
the tubes with an inside diameter of 12.7 mm. Both a 
"direct method," and a "modified Wilson method" were used, 
and are discussed below. 
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Summary of Tubes Tested and Their 
Hea t -Trans fer Performance 
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TABLE II (continued) 
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The Direct Method 



1 . 



For the "direct method," a smooth tube was manufac- 
tured with six wall thermocouples inserted around its perim- 
eter, During data runs under vacuum and atmospheric 
conditions, readings from these six thermocouples were 
obtained and averaged, in order to determine an average tube 
well temperature for each heat flux. It should be noted 
that the wall temperature was guite non-uniform around the 
tube perimeter. Georgiadis [7] reported up to an 18-K drop 
in the wall temperature from the top to the bottom of the 
tube. Evaluation of the uncertainty resulting from this 
temperature variation was considered beyond the scope of the 
investigation. The average wall temperature was then used 
along with the cooling water inlet and outlet temperatures 
to determine the log-mean-temperature difference (LMTD) . 
With values for heat flux and LMTD known, the value of the 
water-side heat- tra nsfer coefficient was obtained. This, in 
turn, was used to determine the leading coefficient for the 
Sieder-Tate- ty pe equation shown below: 



Nu = C Re 0 ' 8 TV 1 ^ 3 



„ 0.14 



l‘w 



+ B 



(5.1) 



A more detailed step-wise solution procedure is given by 
Georgiadis [7]. For the previously-mentioned tubes of 
15.0 mm and 19.05 mm 0. D. , Georgiadis determined the coeffi- 
cient C to be 0.0635. The constant B = 26.4 is an addi- 
tional parameter found to improve the fit of equation (5. 1) . 
Note that C = 0.0635 is about 2.5 times greater than the 
well-known Sieder-Tate coefficient of 0.027 [28] for long 
tubes with smooth inner walls. This is due primarily to the 
insert used to enhance the water-side heat-transfer 
coefficient. 
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2 . 



Th e Modified Wilson Method 

A complete description of this method is provided in 
Appendix B. Briefly, this method uses a Sieder-Tate-t ype 
equation for the water-side heat-transf er coefficient, and 
an equation (discussed below) for the steam-side heat- 
transfer coefficient. Both equations contain leading coef- 
ficients that must be determined iteratively while being 
used in a linear equation to generate the Wilson plot. The 
iterations were performed on the HP 9826A computer mentioned 
previously in Chapter III, using program "WILSDN6." 

The "modified Wilson method" used during this inves- 
tigation consisted of one minor modification from that of 
Georgiadis. For this thesis, a Fujii-type [ 29] equation was 
used instead of the Nusselt-type equation for the steam-side 
coefficient. This modification was made to account for the 
small steam velocity (about 1 m/s) that was pcesent during 
the runs. For this reason, use of the Fujii-type equation 
appears to be more accurate than the use of the Nusselt-type 
equation. This modification resulted in a slightly higher 
value (up to 3 percent) for the leading coefficient C. The 
program "WILS0N6" described earlier, allows an option for 
selecting either the Fujii-type or Nusselt-type equation for 
the steam-side coefficient. 

Using a Nusselt-type equation for the steam-side 
coefficient, Georgiadis [7] found the leading coefficient 
(for equation (5.1)) C to be 0.071, with the B value set 
equal to zero. This C value is about 10 percent higher than 
the value (C = 0.0635) obtained by the "direct method," but 
since the "direct method" is generally felt to be more reli- 
able, the values of 0.0635 and 26.4 were used for the 
constants C and B respectively for this thesis effort with 
tubes of 12.7 mm inside diameter (see Table I). 



3. 



Wate r-S ide C oef ficie nt s For Thin-Wall Tubes 

Due to the thin (i.e., 0.5 mm thickness) tube walls 
for tubes 39 thru 42, it was not possible to manufacture an 
instrumented tube. For this reason, the "modified Wilson 
method," together with the Fujii-type equation for the 
steam-side coefficient, was used to determine the leading 
coefficient for equation (5.1) for both the copper tube 
(tube number 41) , and the stainless steel tube (tube number 
42). For the copper tube, a value of C = 0.0756 was 
obtained, while for the stainless steel tube, a value of 
C = 0.0688 was obtained. Although the inner tube diameters 
for both tubes are the same (i.e., 13.5 mm) , the much higher 
thermal conductivity of copper, which results in circumfer- 
ential heat conduction through the tube wall, leads to its 
higher value for C. A summary of the leading coefficients 
used in equation (5.1) for all the tubes is presented in 
Table III. 

C. F.EPEAT ABILITY OF DATA 

Since the reliability of data taken was of vital impor- 
tance, all data runs were repeated, as Georgiadis [7] did, 
at least three times on different days. The computed steam- 
side coefficients for similar conditions (i.e., same tube 
and the same operating conditions) from different days 
agreed to within ± 5 percent. Additionally, data runs were 
performed on two finned tubes (tubes number 6 and 17), under 
similar conditions to verify the repeatability with data 
taken by Georgiadis. Figure 5.1 shows that the experimental 
steam-side heat-transfer coefficients of Georgiadis and 
those obtained during this investigation agree to within ± 5 
percent. The curves shown in this figure (and all other 
figures of a similar nature that follow) are the least- 
squares-fit curves according to the following equation: 
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n 



(5.2) 



= a AT b 

where a and b are experimentally determine! constants. 
These curves were selected instead of those based on the 
steam-side coefficient versus heat flux since the uncer- 
tainty band on AT is more uniform than the uncertainty band 
based on the steam-side coefficient. A brief discussion of 
the uncertainty bands is presented in Appendix D, along with 
a listing of the uncertainty analysis program "0HA6" and a 
few sample runs. 

D. EFFECTS OF FIN SPACING AND FIN HEIGHT ON PERFORMANCE 

Data were taken on eight tubes (see Table II) with fins 
of rectangular shape (tubes number 26 thru 33) under vacuum 
(approximately 85 mmHg) and atmospheric conditions. Four 
tubes had a fin height of 0.5 mm, and fin spacings of 
1.0, 1.5, 2.0, and 4.0 mm, while the remaining four tubes 
had a fin height of 1.5 mm, and the same sequence of fin 
spacings. Georgiadis [7] reported data on rec-tangular ly- 
shaped fins with the same fin spacings, but with fin heights 
of 1.0 mm, and 2.0 mm. A complete listing of tubes tested 
by Georgiadis is also presented in Table II (tubes 1 thru 
25) . 

Data reductions were performed on the HP 9926A computer, 
using a program named "DRP6 . " This program was essentially 
the same one used by Georgiadis ("DRP4") , with modifications 
to allow for variations in tube diameters, and tube 
material. Data obtained for all tubes tested during this 
thesis effort (tubes number 26 thru 44) are presented in 
Appendix C. 
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1 • Effe ct s of Fin Spac ing 

This section presents results showing the variations 
of the steam-side heat-transfer coefficient with heat flux 
having fin spacing as a parameter. Figures 5.2 and 5.3 
present data for tubes with a fin height of 0.5 mm under 
vacuum and at atmospheric pressure respectively. The smooth 
tube data and the theoretical prediction of Nusselt are are 
included for comparison. In a similar way. Figures 5.4 and 

5.5 present data for tubes with a fin height of 1.5 mm. 
These figures show that the best heat-transfer performance 
occurs for the tubes with a fin spacing of 2.0 mm, while the 
worst performance occurs for the tubes with a fin spacing of 
4.0 mm. Note that these trends are the same both under 
vacuum and at atmospheric conditions. 

The comparison of finned tubes is made through the 
enhancement ratio. This ratio is defined as the steam-side 
heat-transfer coefficient of a finned tube to that of the 
smooth tube (same diameter as the finned tube root diameter) 
at the same heat flux. As was shown by Georgiadis [7], the 
uncertainty in calculating the steam-side heat- transfer 
coefficient increases with decreasing heat flux; so, heat 
flux values of 0.35 and 1.0 MW/m 2 were chosen for vacuum 'and 
atmospheric conditions respectively. Cross plots of the 
enhancement ratio versus fin spacing are shown in Figures 

5.6 and 5.7 for fin heights of 0.5, 1.0, 1.5. and 2.0 mm 
(curves representing e = 1.0 and 2.0 mm are from Georgiadis 
[7]). Enhancements as high as 4.8 and 6.4 can be seen for 
vacuum and atmospheric conditions respectively. Once again, 
fin spacings of approximately 2.0 mm for the tubes with fin 
heights of 0.5 and 1 . 5 mm are shown to provide the optimum 
performance. Data provided by Georgiadis *7] for fin 
heights of 1.0 and 2.0 mm, however, demonstrate an optimum 
heat- transfer performance for fin spacings of 1.5 mm. As a 
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Fin Spacing/ (mm) 

Figure 5.6 Cross Plot Showing Best Fin Spacing Among Tubes 
with e = 0.5, 1.0, 1.5, and 2.0 mm (Vac. Runs). 




CD 



(D 



ID 



r\ 

E 

E 

\ 

cn 

c 

•H 

0 

0 

Q. 

ID 



C 

C\J Ll 






ID 



°3 



OJ 



72 



Figure 5-7 Cross Plot Showing Best Fin Spacing Among Tubes 
with e = 0.5, 1.0, 1.5, and 2.0 mm (Atm. Runs). 



result, there does not appear to be a clear trend for the 
optimum fin spacing as a function of fin height. 

On one hand, the difference in performance between 
the s = 1.5 and 2.0 mm finned tubes, as well as the s = 0. 5 
and 1.0 ram finned tubes is less than the experimental uncer- 
tainty (± 15 percent), so that the optimum spacing may be 
somewhere between 1.5 and 2.0 mm. On the other hand, each 
data point in these figures has been repeated within 5 
percent (i.e., computed steam-side heat-transfer coeffi- 
cient) on different days. Also, as discussed earlier, the 
data taken on tubes 6 and 17 by Georgiadis were repeated 
within 3 percent during this investigation. Therefore, the 
slightly disagreeing trends shown by Georgiadis and by this 
investigator (see Figures 5.6 and 5.7) may be justified by 
the very complex nature of the condensation phenomenon on 
finned tubes. 

The presence of an optimum fin spacing, however, is 
easy to understand. The 1.0-mm fin spacing, which provides 
the largest increase in surface area, performed worse than 
the 2.0-mm fin spacing due to the relatively large conden- 
sate retention angle with the 1.0 mm fin spacing. Visual 
observations showed that the condensate retention angle was 
around 110 degrees for the tube with s = 1.0 mm, while it 
was about 40 degrees for the tube with s = 2.0 mm. The 
additional thermal resistance induced by this layer of 
condensate overpowers the benefit gained from increased 
surface area, so the heat- transf er performance is reduced. 
As fin spacing increases, the condensate retention angle 
decreases, as does the thickness of the condensate film 
between fins on the upper portion of the tube. 3eyond fin 
spacings of 2.0 mm, the decrease in surface area tends to 
overshadow the other effects, resulting in decreased heat- 
transfer performance compared to the tube with the "optimum" 
spacing. 
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Effects of Fin He ight 

As mentioned earlier, fin heights of 0.5 and 1.5 mm 
were examined during this thesis effort, and data for the 
tubes with fin heights of 1.0 and 2.0 mm were available from 
Georgiadis [7]. Table II shows that the general trend is an 
increase of enhancement ratio as fin height increases, which 
results from the increase of surface area provided by the 
higher fins. 

To study the enhancement obtained beyond the 
increase in surface area due to finning, two additional 
columns are provided in Table II. These columns represent 
enhancement ratio/area ratio (Eo/Ar) at both vacuum and 
atmospheric conditions. As can be seen, Eo/Ar values range 
from 0.86 to 2.23 under vacuum conditions, and from 1.43 to 
3.20 at atmospheric pressure. Figures 5.8 and 5.9 present 
plots of of Eo/Ar as a function of fin spacing under vacuum 
conditions, and at atmospheric pressure respectively. The 
relatively small values of Eo/Ar = 0.86 (vacuum) and 
Eo/Ar = 1.43 (atmospheric) for tube number 15 are due to the 
small fin spacing which results in a completely flooded 
tube. On the other hand, the values of Eo/Ar = 2.23 
(vacuum) for tube number 14, and Eo/Ar = 3.20 (atmospheric) 
for tube number 28, show that significant enhancements 
beyond the area ratio are obtainable when the amount of 
condensate flooding is kept small. It appears that surface- 
tension forces are mostly responsible for this enhancement 
beyond the area ratio. 

As discussed in Chapter II, the condensate thinning 
provided by surface-tension effects is related to the pres- 
sure gradient from the fin tip to the fin root. For short 
fins, this pressure gradient is relatively larger than for 
high fins, resulting in improved thinning of the condensate 
film. On the other hand, as fin height increases, surface 
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area increases and the condensate retention angle decreases. 
Both of these effects lead to increased heat- transfer 
performance, but when combined with the poorer thinning 
effect, the Eo/Ar ratio decreases somewhat for the larger 
fin heights. 

E. EFFECT OF FIN GEOMETRY ON PERFORMANCE 

In order to study the effect of fin geometry, data were 
taken on tubes with fins of triangular, trapezoidal, and 
parabolic shapes (tubes number 34, 35, 37, 38), as well as 
on a spirally-finned tube with threads of triangular shape 
(tube number 36). Fin dimensions (i.e. , fin spacing, 
height, and thickness) for these tubes were chosen so that 
the resulting total surface areas would be approximately 
equal to that of the "optimum" tube with rectangularly- 
shaped fins and a fin height of 1.0 mm (tube number 6). In 
this manner, the effects of increased surface area were 
eliminated, allowing a comparison between tubes where fin 
geometry was the only variable. A commercially-available 
"Wolverine" tube (tube number 43), whose fin geometry is 
shown in Figure 3.5(c), was also tested. Dimensions for the 
above fins are given in Table II. 

1 • Effect of Fin Sh ape 

This section examines the effect that fin shape has 
on the heat-tra nsfer performance of finned tubes. The 
performance of tubes number 34, 35, 36, 37, and 43 under 
vacuum conditions is shown in Figure 5.10, while Figure 5.11 
depicts their performance at atmospheric pressure. The 
dashed curve in these figures represents a least-squares fit 
for data taken on the "optimum" tube with rectangularly- 
shaped fins and a fin height of 1.0 mm (tube number 6). In 
both cases, the "Wolverine" tube (tube number 43) gave the 
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q/ CMW/m 2 ) 

Figure 5.11 Variation of Heat-Transfer Coefficient for Tubes 
with Fins of Different Shapes (Atm. Runs). 



worst performance. This was due primarily to the small fin 
spacing (0.7 mm) that resulted in a completely flooded tube. 
Based on the results of Georgiadis [7], the small fin thick- 
ness (0.3 mm) of the Wolverine tube may also be responsible 
for its poor performance. Georgiadis showed that for a fin 
spacing of 1.0 mm the tubes with a fin thickness of 0.5 mm 
(tube number 16) performed 16 percent and 6 percent poorer, 
under vacuum and at atmospheric pressure respectively, than 
the tube with a fin thickness of 1.0 mm (tube number 5) . 
Based on this observation, the small fin thickness is also 
responsible for the poorer performance of the "Wolverine' 1 
tube in comparison, for example, to the tube with fins of 
parabolic shape (tube number 37) , which was also completely 
flooded by condensate. Therefore, the poor performance of 
this tube is probably due more to increased condensate 
flooding and smaller fin thickness than to fin shape. On 
the other hand, the best performance was observed for the 
spirally-finned tube (tube number 36) under vacuum condi- 
tions, and the tube with parabolic fins, where fin spacing 
was 0.7 mm (tube number 37) , at atmospheric pressure. The 
fact that tube number 37 performed well at atmospheric pres- 
sure, but poorer than the spirally-finned tube under vacuum 
conditions was somewhat unexpected. It must be kept in mind 
that these data runs were repeated at least three times on 
different days, thus confirming this trend. One possible 
explanation centers on variable fluid properties, as 
discussed below. 

Under vacuum conditions, the condensate retention 
angle is significantly smaller for tube number 36 (about 30 
degrees) than for tube number 37 (completely flooded) , and 
the saturation temperature is about 50 K lower. Therefore, 
the relatively high heat-transfer coefficient of the 
spirally-finned tube is due mainly to its small condensate 
retention angle. When the condensing temperature is raised. 



80 



both surface tension and viscosity decrease. For example, 
the surface tension decreases by about 10 percent, and the 
viscosity decreases by about 45 percent when the temperature 
is raised from 40 °C to 90 °C (assumed approximate film 
temperatures). Decreased viscosity should result in higher 
heat- transfer coefficients at atmospheric pressure than 
under vacuum conditions for a given tube. Further, it may 
be possible that the effect of viscosity is dependent on the 
fin shape. Unfortunately, no eguations exist to support 
this hypothesis. The effect of surface tension, on the 
other hand, is fairly well understood; the condensate reten- 
tion angle decreases with decreasing surface tension. 
Therefore, the reduction in retention angle (from vacuum to 
atmospheric pressure) for the tube with "parabolic" fins is 
greater than the reduction for the spirally-finned tube, 
since the latter was already small under vacuum conditions. 
Thus, the reduced retention angle on the tube with "para- 
bolic" fins provides an increased heat-transf er performance, 
surpassing the performance of the spirally-finned tube. 

The performance of the remaining two tubes (tubes 
number 34 and 35) is the same under vacuum conditions, but 
at atmospheric pressure, the tube with triangularly-shaped 
fins (tube number 34) outperformed the tube with trapezoidal 
fins (tube number 3 5) by about 10 percent. An argument 
similar to the one above may explain this behavior. 



Effects on Enhancement Ratio 



Referring to Table II, it can be seen that all of 
the above tubes provided significant heat-transf er enhance- 
ments over smooth-tube values. Note that the enhancement 
obtained under vacuum conditions for tube number 36 was the 
same as obtained for the "optimum," rectangularly-shaped 
finned tube with the same fin height (tube number 6) . Under 
atmospheric conditions, however, the tube with fins of 
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parabolic shape (tube number 37) performed significantly 
better than any of the tubes with rectangular fins of the 
same height (1.0 mm) . 

Like the tubes with rectangular fins, most of these 
tubes displayed enhancements above the area ratio, as shown 
by the two Eo/ir columns in Table II. These ratios range 
from 0.94 to 2.95 under vacuum conditions, and from 1.52 to 
4.13 at atmospheric pressure. Like tube number 15, the 
small values of 0.94 and 1.52 for the "Wolverine" tube 
resulted from the large amount of retained condensate, and 
by the very thin fins (about 0.3 mm) as discussed earlier. 
As was the case for rectangular fins, the highest Eo/Ar of 
2.95 and 4.13 were shown for the tube (tube number 36) with 
a relatively small area ratio of 1.37, indicating once again 
the importance of surface-tension effects. 

3 • The £er for ma nee of " Para bolic" F in s 

In order to study the parabolic fin more closely, a 
tube with parabolic fins was manufactured (tube number 38) 
with the same fin height, spacing, and thickness as tube 
number 17 (tested by Georgiadis [7] and this investigator). 
The performance of these two tubes is compared in Figure 
5.12 under vacuum conditions, and in Figure 5. 13 at atmos- 
pheric pressure. Clearly, the parabolic fins provide a 
greater enhancement than the rectangular fins. The reason 
for this, as mentioned earlier, is probably the continuous 
change of radius of curvature (increasing from the trough 
area to fin tip) for the parabolic shape. The resulting 
pressure gradient along the fin height is more favorable 
(than on rectangular fins), which, in turn, leads to 
improved drainage from the fin surface. This effect was 
seen most clearly on tube number 37. Another observation 
(that contributes to improved heat transfer) made on the 
tube with "parabolic" fins is that, although completely 
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Figure 5.12 A Comparison Between Fins of Rectangular Shape 
and Parabolic Shape (Vac. Runs) . 
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flooded, the fin tips (about half the fin height) were 
visible along the majority of the fin circumference (except 
for about the bottom 30 degrees) . 

As discussed in section 2.1, tube number 37 
performed poorer than the spirally finned tube (tube number 
36) under vacuum conditions. The reasons for this are not 
fully understood; so, further tests of these fin geometries 
(i.e., spiral and parabolic) should be performed. 

F. EFFECT OF FIN-METAL THERMAL CONDUCTIVITY ON STEAM-SIDE 

COEFFICIENT 

To investigate the effect of fin-metal thermal conduc- 
tivity on the heat- transfer performance, a stainless steel 
tube (tube number 40) was manufactured with the same fin 
dimensions as the "optimum'’ copper tube (tube number 6) with 
rectangular fins, and a fin height of 1.0 mm. Due to the 

low thermal conductivity of stainless steel, it was neces- 
sary to reduce the tube wall thickness to 0.5 mm to prevent 
an excessive wall thermal resistance. A copper tube (tube 
number 39) with the same dimensions as the stainless steel 
tube was manufactured for comparison. Also, two smooth 

tubes having an outside diameter equal to the root diameter 
of the finned tubes (tubes number 41 and 42) were manufac- 
tured; the first tube was made of copper, and the other was 
made of stainless steel. The results for data runs taken 
under vacuum conditions, and at atmospheric pressure are 
shown by Figures 5.14 and 5.15 respectively. 

As expected, the copper finned tube provided a signifi- 
cant enhancement, while the stainless steel finned tube 
actually reduced the heat- transfer performance slightly in 
comparison to the smooth tube (tube number 42). The thermal 
conductivity of stainless steel is much lower than that of 
copper (i.e., 15 W/m.K compared to 385 K/m.K), so the 
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Figure 5.15 Comparison of Heat-Transfer Performance Between 
Stainless Steel and Copper Tubes (Atm. Runs). 



thermal resistance of the stainless steel fins is signifi- 
cantly higher, resulting in a poorer fin efficiency. This 
is combined with the deleterious effects of retained conden- 
sate (introduced by the use of fins) so that the above 
factors dominate over the beneficial thinning of the conden- 
sate film between fins. The net effect, then, is a slight 
degradation of the heat- transfer performance compared to a 
smooth tube. 

G. EFFECT OF STEAM VELOCITY 

As discussed in Chapter II, most models of film conden- 
sation on finned tubes have neglected the effect of vapor 
shear. For this reason, a decision was made to perform runs 
for the "optimum” tube (tube number 6) with rectangular 
fins, and fin height of 1.0 mm, at nominal steam velocities 
of 2, 4, 6, and 8 m/s. Steam velocity through the test 

section was controlled by the boiler input power, and since 
a steam velocity of 1 m/s was the maximum obtainable at 
atmospheric pressure, the data runs were performed under 
vacuum conditions only. The results of data taken for the 
above-mentioned steam velocities are shown in Figure 5. 16. 

The trend shown in Figure 5.16 indicates that increased 
steam velocity (i.e., increased vapor shear) leads to 

enhancements of up to 10 percent as steam velocity increases 
from 2 to 8 m/s. This is to be expected since the shear 
stress at the liguid-vapor interface tends to pull the 
condensate toward the bottom of the tube, thereby thinning 
the condensate film at the top part of the tube. 

The above-mentioned data were compared to data taken by 
Yau et al. [1] for steam velocities of 0.5, 0.7, and 

1.1 m/s at atmospheric pressure. They reported enhancements 
of up to 40 percent for the steam velocity of 1.1 m/s 
compared to the 0.5 m/s velocity. The rather large 
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difference in enhancements reported by Yau et al. compared 
to those of this investigation was not expected. It should 
be kept in mind, however, that the data of Yau et al. were 
collected at atmospheric pressure (compared to vacuum condi- 
tions for this investigation), where the increased satura- 
tion pressure leads to a lower condensate viscosity (about 
45 percent as discussed in section E. 1) . This, in turn, 
allows condensate to flow more freely so that the effect of 
shear stresses acting on the liguid-vapor interface would be 
more pronounced. Another factor to be considered is the 
slightly larger fin spacing used by Yau et al. (2.0 mm 
compared to 1.5 mm for this investigation). If the fin 
spacing is very small, the vapor shear may not be "felt" on 
the sides of the fins or on the fin bases. This is due to 
the fact that the steam flows around the tube, thus leaving 
nearly stagnant steam between fins. As the fin spacing 
increases though, the fin sides and bases will be more 
"accessible" to steam, therefore, giving rise to the benefi- 
cial effects of vapor shear. Thus, the larger enhancements 
of Yau et al. [1] compared to this investigation may be 
explained by the two phenomenon just described. 

An attempt was made to correlate the data using a 
rujii-type [29] equation given by: 
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(5.3) 



where 



Mu = h af n o 
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(5.5) 



Re 
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= V s p o 
Pf 



F = S D o g f h ffi 

V s z k f AT 



(5.6) 



V = steam velocity 

S u 

h af = steam-side heat-tra 
on actual surface a 

Re fcp = two-phase Reynolds 
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Figure 5. 17 Comparison of Experimental Data to 
Fujii-Type Equation (Equation 5.3). 



Note that the slope of the line predicted by equation 
(5.3) for the experimental data does not match the slope of 
the Fujii equation mentioned earlier, but instead, nearly 
matches the slope of the Nusselt equation for smooth tubes 
given by: 



h = 0.723 



3 2 

k f Pf 3 h F( , 

Uf 0 o AT 



1/4 



(5.7) 



which can be rewritten as: 



Nu 



1<e tp 



1/4 

0.72R F 



(5.8) 



In this manner, the effect of vapor shear, at least on the 
tube tested (tube number 6) , is very small compared to 
smooth tubes [7]. In fact, using the experimentally-found 
d-value in equation (5.3), the steam-side coefficient and 
steam velocity are related as follows: 



Once again, this observation shows that the effect of steam 
velocity on the steam-side hea t- transfer coefficient is very 
small. 

A vertical shift of the least-squares-f it compared to 
the Nusselt theory results from the insreased c-value (1.302 
compared to 0. 728) of equation (5.8). In fact, the ratio of 
the experimental c-value to the coefficient in the Nusselt 
equation (i.e. , 1.302/0. 728 = 1.788) represents the heat- 

transfer enhancement in excess of the area increase due to 
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finning. Therefore, the c-value depends fin geometry, and 
surface tension of the condensate etc.. 

In view of the above discussion, further attempts to 
modify the Nusselt equation to include vapor shear, surface- 
tension effects, and fin geometry etc. may lead to an accu- 
rate correlation for the steam-side heat- transfer 
coefficient for finned tubes. 

H. EFFECT OF INTERNAL AND EXTERNAL ENHANCEMENTS ON THE 

OVERALL COEFFICIENT 

Dp to this point, all discussions have delt with 
enhancements to the steam-side heat-transfer coefficient. 
Even with large steam-side enhancements, however, a large 
thermal resistance through the tube wall, or on the water 
side may dominate, and completely negate the steam-side 
enhancement. For this reason, a decision was made to test 
the "optimum” tube with rectangular fins (fin height of 

I. 0 mm - tube number 6), and the corresponding smooth tube 
(tube number 1) , with and without the insert. The results 
of these data runs, taken under vacuum conditions, are shown 
in Figure 5.18. At a water velocity of about 3 m/s, use of 
the insert in the smooth tube resulted in a 40 percent 
increase of the overall coefficient, as did the use of 
external fins without the insert. By combining the use of 
fins and the insert, an enhancement of nearly 155 percent 
was attained. It must be noted, however, that the insert 
used during these experiments is not suited for use in 
condensers due to the very high pressure drop it creates. 
Instead, Webb et al. [30] showed that the use of internal 
flutes, fins, or circumf erential ribs may be considered to 
increase the water-side coefficient by a factor of up to 2.0 
(compared to 2. 5 obtained by the insert used during this 
investigation). As a result, it appears that the use of 
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external fins togather with internal enhancements, could 
increase the overall heat- transfer coefficient by approxi- 
mately 100 percent. 

I. DEVELOPMENT OF AH EXPERIMENTAL CORRELATION 

Using the 1985 model of Webb et al. [25] (discussed in 
Chapter II) as a starting point, attempts were made to 
correlate the experimental data obtained for tubes with 
rectangularly- shaped fins. Time limitations prevented an 
in-depth effort, but a few simple modifications were made. 
Before discussing these modifications, the general scheme of 
the Webb et al. model will be discussed. A more detailed 
step-by-step discussion is provided in Appendix A. 

First, a value for the parameter £ must be determined 
for the fin geometry in question. With this value known, 
the heat-transfer coefficient for the unflooded fin surfaces 
can be obtained using equation (2.30). Values for the heat- 
transfer coefficients of the unflooded tube area between 
fins, and for the flooded region are determined next using 
equations (2.32) and (2.33) respectively. The above- 
mentioned heat- transfer coefficients are then substituted 
into equation (2.31) to obtain the average heat- transfer 
coefficient for the finned tube in question. 

A careful review of the Webb et al. model revealed two 
deficiencies worth noting. First, in their two-dimensional 
analysis of the flooded tube region, the heat-transfer coef- 
ficient for the flooded root area between fins was applied 
to the fin tips. Figure 2.4 shows that a constant, rela- 
tively thick condensate film is assumed which, when applied 
to the fin tips, results in an underprediction of the heat 
transfer. Use of the Nusselt equation for smooth tubes 
would be more appropriate for fin tips in this region. 
Additionally, the effects of convection in this region are 
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Figure 5.18 Overall Heat-Transfer Coefficient for 
a Smooth Tube and Finned Tube with 
and without Insert (Vac. Runs) . 



to include these effects would make a 



ignored; however, 
solution extremely complex. 

The second deficiency, probably valid only for 
rectangularly-shaped fins, involves the assumed condensate 
film shape between fins on the unflooded tube surface. 
Figure 2.4 shows a depression in the film next to the fin 
base, with an increasing condensate thickness away from it. 
Condensate film shapes observed during this investigation 
continued to decrease away from the fin base as shown by 
Figure 2.1, probably due to the thinning effect of surface- 
tension forces. Based on the assumption of a relatively 
thick condensate film, predictions using the Webb et al. 
model will underpredict heat transfer from this part of the 
tube. 

To use the Webb et al. model for rectangular fins, a 
value of -0.9 was chosen for £ . A value of -1.0, which 
corresponds to a straight line (i.e., straight side of the 
rectangular fin) would seem more appropriate; however, this 
value would lead to zero heat transfer from the fin sides 
using equation (A.1). Webb et al. used values around -0.7 
to -0.8 for their slightly trapezoidal fins; so the value of 
-0.9 was assumed to be reasonable. Using this value along 
with a <j> value of 5.0 (used by Webb et al.) resulted in 
significant underprediction of the experimental data 
obtained by this investigator and by Georgiadis [7], as 
shown in Figure 5.19. 

In an attempt to improve the predictions, values for the 
<J> ratios were calculated using the paralled heat flow paths 
proposed by Owen et al. [13] discussed earlier in 
Chapter II. These <j> values are dependent on fin spacing, 
and should decrease with increasing fin spacing. This trend 
was observed, and for the set of tubes with both fin height 
and fin thickness of 1.0 mm (s = 0.5, 1.0, 1.5, 2.0, and 
4.0 mm - tubes number 4 to 8), resulted in 4> values that 
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ranged from 1.7 (s = 0.5 mm) to 1.2 (s = 4.0 mm) . Using the 
Febb et al. model once again, but with a constant average 
value of 1.5, gave the results shown in Figure 5.20 for the 
above-mentioned tubes (tubes number 4 to 8) . This Figure 
also shows a plot of the Webb et al. model using the vari- 
able (i.e., 1.7 to 1.2) which made almost no difference. 

Once again, the experimental data, also shown in Figure 5.20 
were significantly underpredicted. 

In order to correct for the assumed thick condensate 
layer between fins, the work of Fujii et al. [31] was used. 
Their paper delt with smooth tube enhancements from coiling 
a wire around the tube circumf erence , and provided an equa- 
tion to account for the thinning effect between wire wraps 
as shown below: 



Nu f 

Wu 8 
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(s + 1 1 



4 (l + A) ‘ 
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1/4 



(5. 10) 



where 
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2 a 
Pf R 
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(5. 11) 



(5. 12) 



D w = wire diameter (taken as the fin height 
for use in the Webb et al. model) 



Adding this correction factor to the heat- transfer coeffi- 
cient for the unflooded root area between fins made a 
significant improvement for the group of tubes mentioned 
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earlier (tubes number 4 to 8) as shown by Figure 5.20. 
Nonetheless, experimental data were still underpredicted. 

In summary, it should be kept in mind that this model 
was not developed for rectangular fins. The most signifi- 
cant discrepancy, however, appears to be on the prediction 
of heat transfer through the flooded area. This model 
predicts only minimal heat transfer through the flooded 
portion of the tube, whereas Fanniarachchi et ai. [2] 
reported enhancements at least equal to the area ratio for 
completely flooded tubes. Additionally, there may be other 
phenomena which have not been taken into account (such as 
intense convection in the condensate} which, if included, 
may improve the predictions. Nonetheless, this model shows 
promise and should be pursued further. 
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VI. CONCLUSIONS AND RECOMME NDATIO NS 
A. CONCLUSIONS 

1. The use of fins may lead to significant enhancements 
of the steam-side heat-transfer coefficient. Under 
vacuum conditions, enhancements as high as 4.0 (tube 
number 36) were realized, while enhancements of up to 
6.9 (tube number 37) were realized at atmospheric 
pressure. These enhancements, which were up to three 
times greater than the area ratio (finned tube 
area / smooth tube area) are due to the thinning 
effect of surface-tension forces. 

2. For rectangularly-shaped fins with fin heights of 0.5 
and 1.5 mm, the 2.0 mm fin spacing appears to be the 
optimum among the fin spacings tested (1.0, 1.5, 2.0, 
and 4.0 mm) . This was the case for both vacuum and 
atmospheric conditions. For similar tubes with fin 
heights of 1.0 and 2.0 mm, Georgiadis *7] found the 
optimum spacing to be 1.5 mm. 

3. The use of fins on a stainless steel (low-thermal- 
conductivity- metal) tube degraded the heat-transfer 
performance in comparison to its corresponding smooth 
tube. In this case, the poor fin effeciency combined 
with the deleterious effects of retained condensate 
dominated over the beneficial thinning of the conden- 
sate film between fins, resulting in a net reduction 
of the heat- transfer performance. 

4. Fin geometry affects the extent to which the surface- 
tension forces thin the condensate film on the fin 
surfaces. Parabolically-shaped fins were shown to 
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maximize the thinning effect of surface-tension 
forces. 

5. Increasing steam velocity (i.e., vapor shear) from 2 

to 8 m/s increased the steam-side coefficient by 
about 10 percent (under vacuum conditions) , in 

contrast to Yau et al. [1]# who showed a 40 percent 
increase when steam velocity was increased from 0.5 
to 1.1 m/s (at atmospheric pressure). 

6. For the "optimum" tube with rectangular fins and fin 

height of 1.0 mm (tube number 6), the overall heat- 
transfer coefficient without the insert showed an 
enhancement of 40 percent over the corresponding 
smooth tube. With the insert in place, the finned 

tube showed overall enhancements of about 155 

percent. 

7. Host theoretical models developed thus far to predict 
the heat- transfer performance of steam condensation 
on horizontal finned tubes neglect one or more 
phenomena of major importance. Thus, there are no 
models currently available to accurately predict the 
performance of finned tubes used to condense steam. 

B. RECOMMENDATIONS 

1. Make further refinements to the Webb et al. [25] 
model in an effort to predict the experimental data 
within a reasonable accuracy. 

2. Manufacture and test additional finned tubes with 
parabolic and spiral fin shapes with different dimen- 
sions (i.e., fin spacing, height, and thickness). In 
this manner, determine the optimum dimensions for 
each fin shape. 
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3. Manufacture and test a series of tubes with varying 
tube diameter to study the effect of diameter on the 
heat-transfer performance. For this purpose, select 
the optimum, rectangular-shaped fins (s = 1.5 mm, 
e = 1.0 mm, and t = 1.0 mm). 

4. Collect data on finned tubes fitted with drainage 
strips made out of both solid and porous metal to 
determine the enhancements they provide. 

5. Continue tests on finned tubes with tube-metal 
thermal conductivities between those of copper and 
stainless steel, such as copper-nickel alloys and 
aluminum. 

6. Conduct data runs at atmospheric pressure with steam 

velocities of 0.5, 0.7, and 1.1 m/s, using tube 

number 6 to compare with the data of Yau et al. [1]. 

7. Perform data runs using different fluids with widely 
varying fluid properties (i.e., surface tension, 
viscosity, etc.) to further study these effects, and 
to compare with existing theoretical models and 
experimental correlations. 
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APPENDIX A 



PROCEDURE FOR OSE OF THE WEBB ET AL. [25] MODEL 

As proposed by Webb et al., solving for the steam-side 
heat-trar.sfer coefficient proceedes as follows: 

1. First, a value for £ must be determined. This 
parameter, as shown by figure 2.4, is a measure of 
the fin aspect ratio (e/t) , and it must be deter- 
mined using an iterative method. It is a function 
of fin geometry only, and should match the fin side 
profile as closely as possible. Webb et al. [25] 
provide details of this method, and no further 
discussion is provided on this since a true ^-value 
does not exist, as defined by Adamek [27], for a 
rectangular-shaped fin. 

2. With the value of known, determine the heat- 
transfer coefficient for the fin sides in the 
unflooded portion of the tube. This value is 
obtained using the following equation "25, 27]: 



h, = 2.140 

r s 



o h f p 0 m S m pf (£+ 1) 

U f kf AT 7J~+ T) J 



1/4 



(A.1) 



3. Determine the heat-transfer coefficient for the 
unflooded root area between fins using the Nusselt 
equation written in terms of the film Reynolds 
number. For this purpose, use the following itera- 
tive method: 
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Calculate the condensation rate (m) using: 



m = n A f( . AT 

hfg. 



(A. 2) 



Assume: 



m = m 



(A. 3) 



Calculate the film Reynolds number by: 



Re - 4 

Uf r r s + t h ) - t h 1 



(A. 4) 



Calculate h r using: 



h = 1.514 
r 



2 ^- 1/3 

Re , 



kf 3 Pf z g f 



(A. 5) 



Calculate the base heat flow by: 



°r " \ V AT 



(A. 6) 



Solve for m using: 



m = m + _I 
r h 



fg 



(A. 7) 



Repeat steps c thru f until convergence is 
obtained . 



Determination of the heat- transfer coefficient in 
the flooded part of the tube requires that a 
computer solution be used to solve the two- 

dimensional conduction problem of heat flow through 
the fins and retained condensate. This solution 
provides a ratio (<t>) of actual heat flow into the 
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The fin efficiency is determined next using the 
following equation: 
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6. The heat-tra nsf er coefficient for the entire tube 
then, is obtained by substituting the above guanti 
ties into the following eguation: 
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ow 
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n 

o 




^1+nh. 

A sf f 



A ft 



A sf 
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TT h TT 



(A. 11) 
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APPENDIX B 

MODIFIED WILSON METHOD 



The "modified Wilson method" calculates the leading 
coefficient for the Sieder-Tate equation. This method 
assumes a form of correlation for both the water-side and 
steam-side heat-transfer coefficients, each of which 
contains a coefficient to be determined by iteration. In 
the past, a Nusselt-type equation was used for the steam- 
side heat-transfer coefficient as shown below '7]: 

3 1/3 

k f Pf ^ P f ” P v ^ h f & 

S7^\T5 J < B * 1 > 




Equation (B. 1) results in a 6 value of 0. 655 for conditions 
where no vapor shear is present. The presence of vapor 
shear generally results in a higher value, which must be 
determined iteratively. The use of this equation should, 
therefore, be restricted to relatively low vapor shear 
conditions (i. e. , vapor velocity less than 1.0 m/s). Note 
that a Nusselt-type equation does not represent the "correct 
form," as it contains no terms involving vapor shear. To 
avoid this deficiency, a correlation developed by Fujii and 
Honda [29] which takes the effect of vapor shear into 
account was used, and is shown below: 
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Rewriting equation ( 3 . 2 ) to express h as a function of heat 
flux and vapor velocity results in: 







i/4 


= 6 


R h fR 











3 3 3 

"f D o Pf V s 



„ 1/8 
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(B.3) 



The constant g in eguation (B.3) is to be determined by 
iteration. A Sieder-Tate- type equation is used for the 
water-side hea t- transfer coefficient, as shown below: 



h, n, „ 0.8 „ 1/3 

_i L = c Re Pr 

kf 



J w 



o. m 



= c n 



(B.4) 



Equations (B.3) and (B.4) are then substituted into the 
equation for the overall heat-transfer resistance given 
below : 



1 + + £* 



"c A o N A 1 h 0 A n A n 



(B.5) 



A linear eguation used to generate the Wilson plot is 
obtained, and shown below: 



Y = 



* + i 
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(B. 6) 



v = r 
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(B . 8 ) 
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The parameters X and Y are determined from fluid-property 
values and heat-flux measurements taken during data runs 
using an uninstrumented smooth tube. Iterations between the 
coefficient C, and the coefficient B , are continued until 
convergence of the coefficients (between two successive 
iterations) within 0.1 percent occurs. The slope of the 
Wilson plot generated is the reciprocal of the desired 
Sieder-Tate coefficient. 



APPENDIX C 



LISTING OF RAH DATA 

The following pages contain raw data obtained for tubes 
number 26 thru 44 under vacuum conditions and at atmospheric 
pressure. Raw data for tubes number 1 thru 25 are presented 
by Georgiadis [ 7 ]. 
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Tube Number : 3 2 T ub e N i un b e r : 3 3 

File Nume ; F82k- , .:6 File (JcMne ; F32N68 
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Tube Number: 67 lj.be Number: 27 

File Name: F 37 P‘3 7 File Name: F37A103 

Pressure Condition: Vacuum Pressure Condition: Atmospheric 
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T ube Number ; j8 1 ube Numoer: :!8 
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&PPEN2II D 

UNCERTAINTY ANALYSIS 



There is always an uncertainty associated with any meas- 
urement which is dependent on the measuring-device accuracy, 
as well as on the operator’s experience. Numerical data 
collected during this thesis effort were used together with 
theoretical formulations, so final values of the steam-side 
heat- transfer coefficient may be distorted due to uncer- 
tainty propagation during calculations. In cases where the 
final results show large uncertainties, it may be unwise to 
draw any conclusions. Instead, the apparatus and/or the 
measuring techniques may have to be modified. 

Estimates of uncertainties for this investigation were 
obtained using program "UNAS" which is listed at the end of 
this appendix, along with selected uncertainty evaluations. 
This program determines the uncertainties using an equation 
proposed by Kline and McClintok [32] shown below: 




where 

W R is the uncertainty of the desired dependent 
variable 

x-l , x 2 , ...., x n are the measured (independent) 
variables 

, w 0 , ...., w n are the uncertainties in the measured 
variables 
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A complete discussion covering the de velopement of 
uncertainty analysis used for this investigation is given 
Georgiadis [7]. 



th e 
b y 
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M90 ENTER 'jP : l e : I f g . Inn 
1 ! 3 T IF It *3 = U THEN ENTER ?File:Dd 

1192 IF Ifg-i THEN ENTER *Fi le : Da .Dd. Dd j 

1200 FOR I*! TO las 

1 2 M3 ENTER s»F * ie : Bvo i , Samp . Vtran.£tp.E<M . r m . 7c i .Tcn.Phg. P water 
1220 NEXT I 
*230 Enf*E<0> 

1240 I c P re- 1 THEM 
1250 BEEP 

*260 3 RINT USING " 1 OX , " M °ressure Condition: Vacuum <11 if? a) 

1270 ElSE 

* 230 PRINT USING " MIX , ""Pressure Condition. Atnospneric M01 irPa> 

1230 END IF 
1300 PRINTER IS : 

1 ?; 0 BEEP 

1320 C 'RINT USING "4x . ""Select xudc* .nail r yoe: G-truc* l-tmn 

i 330 INPUT It: 

? J4|J IF Itt = u THEN 
1 9 50 Do a .01905 
1360 D i m Ai • 27 
1270 ELSE 
1 280 Do-. 01 45 
I 29u Di- .0135 
1400 END IF 
1410 PRINTER IS 1 
1 420 BEEP 

i 43U PRINT USING "4X . ""Se lec t material code: 

1440 PR 1 14 T USING "4X . ""G-Cll, • -SS 

1450 INPUT Inc 
1460 IF lmc-0 THEM 
14 70 k'c-385 

1 47 1 Dkc-10 

1480 Ci-. 07’ 

1 490 Ac-0 
1500 ELSE 

1 5 1 0 Kc-15 

1511 Dice- 1 
1520 Ci -.0688 
1530 Ac=H 
•540 EMD IF 
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1550 P R INTE P [3 1 0 I 
1 560 Ts 3 FN Tvsv < Enf > 

1570 PRINT USING '* 1 GX . " M S team Tenoerature * "" . 3D . 2D . *"* <Deg C> : 

I 3 

!53D p RIN r USING '* * OX , ‘‘‘"Water Fiou Pace <7) * ’’W 30 . 2D" : Fn 

1 5*30 Dtci-.Qn: 

1 600 Dtco 3 .005 

1610 BEEP 

I 620 Demf - I . OE -6 

1 630 D t s 3 SGR < ( < C < I 3 ) •Enf C *4 *C ( 4 ) •Ent * 3 ) • Den f ) 2) 

1640 T 3 < Tc i + Tco > /2 ! FILM TEMPERATURE 

1650! UNCERTAINTY [N THE COOLING WATER 
1660 Drho - .5 ! ERROR IN WATER DENSITY 
1670 Dmt • . 009 ! ERROR IN NASS F-.QW KATE 
1630 Rho-FNRho < T > ! WATER DENS ITT 

1690 Mf-I .')48USE-2+6.30332E->Ffn ! HASS FLOW PATE OF COOLING WATER 

1 700 ! CORPECT MF FOP THE TEMPERATURE EFFEO T 

1710 rtf-Mf-M .0365-: ^6S4 «E->t C i-5.252E-6- t ci 2>/ 3954 1a 

•720 Ai 3 (PI-Dl 2>/4 ! TUBE INSIDE CROSS SECTION AREA 

1730 Ddi 3 . 000025 

1740 'Dai J <PI-0di 2>/4 ! ERROR OF INSIDE TUBE CROSS AREA 
1750! COMPUTE THE WATER VELOCITY 

1 760 Vw*Mf / < Rho-A i ) * WATER VELOCITY , 

1 770 PRINT USING " I OX . "‘’Water Velocity - M ‘\Z.QD." M <m/s> ;V 

WOO! CORRECT OUTLET hATER ~E-‘ 4 P . r CR THE* FIXING CHAMBER EFFECT 
1 730 IF I nn 3 1 OR Inn-5 ' P EM 'co 3 Too- . Qi)*-Uiu 2 
1600 IF lpn*0 THEN " co s "co - 1 - .00 I IS* .001 * % hj 2“) 

1810 TmTciWco>*,5 • FILM TEMPERATURE 

1320! COMPUTE THE ERROR IN *ATER VE-JCITY 

1330 Dvu-V.u*SQR< ( Dmt /Mf > 2 + ( Drno/Rho > 2+<Dai/Ai> 2> 

iSaOi UNCERTAINTY IN THE REYNOLDS NUMBER 

1350 Mw-r NMu< T ) ! WATER VISCOSITY 

1 360 Dmu“6.E-b ! ERROR OF WATER VISCOSITY 

1370 Re* < Rno*Vij*D i ) /Mu 

1330 Dr e*Pe*SQR < < Q r no/Rho > ‘2 f < D'>u/ Vu ) 2+<Ddi/Di) 2+<Drou'Mu) 2> 

1390! UNCERTAINTY IN THE HEAT TRANSFERRED 
1300 Cp*u 3 rNCpu< T) 

1910 Q=Mf • < Tco- Tc ; > *Cpw 
1920 Dcpu-8 

J 330 Dq-Q»SGR< < Dmt /Mf ) m 2 + < < D tco/ < Tco- To i ) > > 2 *< < 0 tc i/ < Tco- T c i > > > 2 + < Dcpu>/Cpu J 2 
) 

1340! UNCERTAINTY IN r HE -»EAT FLUX 
1950 D l =.0005 ! ERROR IN TtjBE LEUGKT 
1360 Ddo«. 0*10025 

1370 L 3 , 13335 ! CONDENSING T URE LENGTH 

1930 Qd 3 Q/(PI-Do-L) ! HEAT FLUX 

1990 PRINT USING " T OX . ""Heat -lax 3 . Z . 3DE . "" <W/n 2) : 

Qp 

2000 Dap 3 ( Gp/P I > ♦SQR < ( D q/ 0 ) 2*< Dao/Do ) "2 + < D 1 /L > 2) 

2010 Lm td 3 < Tco-Tci)/L0G( ( Ts - T»; i > / < 7s- Tco ) ) 

2020 Uo-Op/Lm td ! QVERALl. hEAT TRANSFER COEF . 

2030 A1 3 Dts 3 < Tc i - Tco > / < ( 7s-Tci>*< Ts-Tco) 3 L0G ( ( Ts-Tci)/*Ts-Tco> » > 

2040 A2 3 Dtci/< (Ts-Tci > -lOG* < Tj-Tci )/< Ts-Tco > > ) 

2050 A3 3 D tco/ < < Ts-Tco ) -LOG < ( Ts- Tc i > / < Ts-Tco) > > 

2060 D Imtd’Lm td*SQR< A I '2 f A2 '2*93 '2 ) 

2070 Du.o“Uo*SQR ( < Dap/lip ) 2*< Dimtd/Lmtd) *2> 

2030 M=Mu 

2090 1 1 3 < T *27? . i 5 > /273 . ! 5 

2100 Kw*FNKu< I 1 ) 

21 iO Ac-0 . 1 INTERSCEPT PROM SIECcR PROGRAM 
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2120 
2 no 

2140 
2 s 5.0 
2160 
2 170 
2180 ! 
2190 
2200 
2210 
2220 
2230 
2 240 
2250 
2260 
2270 
2280 
2290 
2200 
2310 
2320 
2330 
2 340 
2350 
2760 
2370 
2280 
2330 
2400 
24 i o 
2^20 

2 <* 30 ! 
2440 
2450 
2460 

1 > > > 
2470 
2490 
2490 
2500 
2510 
2520 
2530 
? 5<*0 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 

2 bviO 

2640 
2650 
2660 
2670 
2680 
2 £ 90 
2700 



LI - . 060325 ! LENGTH OF IJMFINNE? LEFT PART OF TUBE 

L2* . G84925 ! LENGTH OF jNFINNED RIGHT PAK T OF TUBE 

01 *.01 905 

D2*. 015875 

Pr*Cpu*Mw/> Cw 

Muw-FNMuu < T ) 

UNCERTAINTY OF INSIDE HEAT-TRANSFER CQEFF . 

C*-1. 

Hi=<Ku/Di)*<Ci*Ke'.S*Pr . 333 ♦Of -Ac ) 

Dt i *Q/ <PI*D i*<L+L 1 *Fel +L2*Fe2 >*Hi ) 

Cf c = < Muw/FNMuw i T + Dt ; ) > '.14 
IF AB$< (Cfc-Cf >/CfO> .01 THEN 
Cf - < Cf +Cf c )* . 5 
GOTO 2200 
END IF 

PI *PI *< D i + D1 ) 

BWD1-bi )-PI-(Di^D! ) * .5 
M 1 - ( H i *P1 / < kc*3l ) ) .5 
P2*PI*<Di-D2> 

B2* <D2 -D i >*PI * < D i +D2 > * . 5 
M2*< H i -P 2/ ( <c-B2 ) > ,5 
F e 1 - F N 7 an h < N 1 *L 1 > / < M 1 *L 1 > 

Fe2*^N7anh ( M2*L2 > / ( M2 *L2 ) i 

D t c *Q / f P I * D i * < L + L *' * F *? 1 ♦ L 2 * r e 2 > • H i ) 

IF ABS( (Dtc-Dti )/Dtc)>. 01 THEN 2200 

DImj-.COIO ! ERROR IN WATER THERMAL CONDUCTIVITY 

Dc i - . 0005 ! ERROR IN SIEDER-TATE COEFFICIENT 

Dor- . 05 ! ERROR IN PRAND7. NUMBER 

Dct *8 . E-6 

A4= , i4*Qcf /Cf 

Dru-hi-SGR( (Dku/Ku) 2*<Ddi/Di > 2+< .3*Dre/Re> 2+< . 333*0pr/Pr ) 2 f <Dc:/Ci) 2- 

UNCERTAINTY OF OUTSIDE HEAT- TRANSFER CHEF-. 

Rw = Do-lOG < Do/D i ) / ( 2*kc ) ' WALL RESISTANCE 

Ho= 1 / < ( !/Uo)-<Do-i_/<Di*<L' , 'Li-rel + L2-Fe2)-Hi ) ) -Ru ) 

Dru*Ru*SQR( ( Dao/Do > ‘ 2 + < Df c / >< c > " 2 + ( Ddo/ < Do-LOG < Do/D i ) > > 2 * < Dd i / ( D i *lOG ( Do/D 
2) 

A5* 1 /'Jo-Ru- ( Do/ ( D i *H i ; ) 

A6-Duo/ ( Uo' 2*A5 > 

A 7 -Oru/AS 

A8*< ( Do/ < D i *H i ) )-(D hi/Hi > >/A5 
PRINT 

Dno = Ho*SQR < A6 2+A7~2*A8 2> 

CALCULATE THE 7. UNCERTAINTY IN Ho 
Pr ho*Dho* ( 00 /Ho 

CALCULATE THE 7. UNCERTAINTY IN REYNOLDS NUMBER 
Prre*Dre*1 00/Re 

CALCULATE THE 7. UNCERTAINTY IN MASS FLOW RATE 
Prmt*Dm»*lOO/Mf 

CALCULATE THE 7 UNCERTAINTY IN HEAT TRANSFER 
Prqp=Dap* 1 00/Qp 

CALCULATE THE 7. UNCERTAINTY IN LMTD 

Primtd*Dimtd*lOO/Lmt<j 

CALCULATE THE /. UNCERTAINTY IN Ry 

Prry=Dru*100/Ry 

CALCULATE THE X UNCERTAINTY IN OVERALL HEAT TRANSFER CDEF . 

P rue* Duo* 1 00 /Uo 

CALCULATE THE 7. UNCERTAINTY IN INSIDE HEAT TRANSFER COE^F . 

Prni=0hi*100/H: 

PRINT 

PRINT USING M T0X. ""UNCERTAINTY ANALYSIS: 
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2710 PRINT 



2 7 n 


PRINT 


USING 


M 1 OX ( 


VARIABLE 






P ERCENT 


UNCERTAINTY 


2712 


PRINT 
















2 720 


PRINT 


USING 


" 1 ox , 


, M *’ M ass Plow Ra t e . 


rid 






?D . 2D . “ : ? rnf 


2730 


PRINT 


USING 


•*iox, 


, ““Reynolds Nmnoer 


. Re 




" •* 


3D . 2D . “ :P r r e 


2 7 40 


PRINT 


USING 


“ i ox , 


, “"Heat Flux, q 








. 2 . 2D . “ : Prop 


2750 


PRINT 


USING 


"10X. 


. ““log-Hean- Ten Di 


f f . 


LtfTD 




. 1 . 2D“ : Pr In td 


2 7S0 


D RINT 


USING 


" 1 OX , 


/'“Wall Resistance 


, Ru 




tl •< 


’ . DD . 2D . “ : P r ru 


2770 


PRINT 


USING 


M 10X , 


. '"’Uveral 1 H . T . C . . 


Uo 




M • < 


’ • DD . 2D . M : Pruo 


2730 


PRINT 


USING 


•• i ox , 


, ““Water- 3 i ae H . T . 


C. . 


Hi 


• •II 


3D .2D . “ :Prhi 


2730 


PRINT 


USING 


•*iox 


. ““S team-S i ae H.T. 


C. , 


Ho 


*• #f 


, 3D . 20 . “ ; Prno 



2300 END 

2810 DEF FNMuw<T> 

2820 A*24 7 . 8/ < T ♦ 1 33 . 1 5 > 

2830 Muu=2.4E-5-lO~A 

2840 RETURN Muu 
2850 FNEND 
2380 DEF r NTann( X 1 
2870 P-EVP<X> 

2880 Q-EXP< -X) 

2890 Tann* < P-O) / ( P-^-'j > 

2300 RETURN Tann 
2910 FNEND 
2920 DEF FiNKuf T 1 ) 

2930 Ku*- .92247 + T I * < 2 . 3395- T ’ * < 1 . 800 7-T I -< .52577- . 07344* T1 > > > 

2940 RETURN ku 
2950 PNEND 
2980 DEF r NMu<T> 

2970 P-247 .3/< T*I33. !5> 

2980 Mu-2 . 4E-5- 1 0 8 
2990 RETURN Mu 
3000 FNEND 
3C. 1 0 DEF FNRho' T) 

2020 Rho-999 . 52948 + T * < . 0 1259- 1 * < 5 . 4825 1 3E- 3- T • 1 .234147E-5) > 

3030 RE 'URN R no 
2u40 FNEND 
3050 DEF FNCdu<T> 

3080 Cou* <4.21' 20858- T •< 2 .26326E- ’- T-< 4 . 4236 12-5*2 . 7 1 429E-7»T ) ) > » 1 000 

3070 RETURN Cpu 

3080 FNEND 

3050 DE” rNTvsw^Em-i > 

3100 COM /Cc/ C< 7 > 

3110 T-C<C> 

3120 FOR I »i TCI ? 

3130 T-T + C < I> *Emf I 
J I 40 NEXT : 

3150 RETURN T 
3180 FNEND 
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DA I A FOR the uncertainty 


ANAL YS I 


C i 


File Name: 




F32V26 


Pressure Condition: 




Vacuum < ’ 1 Kr‘a.> 


Stean Temperature 


= 


a 8. A 7 ( Deq O 


Hater F iow Rate < '/< > 


= 


20.00 


Hater Velocity 


- 


i . ! 6 in/s) 


Neat Fiux 




2.111E + 05 (W/ir. 2) 


UNCERTAINTY ANALYSIS : 






CAR TABLE 




PERCENT UMCEFTAIN 


Naas Flow Rate. Md 




5 . at 


Reynolos Number, Re 




5 , A 8 


Heat F i u.x . q 




i . 7 A 


i - . _m _ - , _ r‘i : r 1 Min 

l_ -» i 'Can t <r l.' i i T * i-. . t ♦ 




0 o? 

v • U fc. 


Nall Resistance. Ru 




2.6? 


Overall H . T . C . . Uo 




1 .75 


Water -Side H.T.C.. Hi 




4 . 46 


S t earn - S : de H . T . C . . Ho 




72.53 



;uv 



h r OR THE UNCERTAINTY ANALYSIS 



- i i e 


■ iane : 


"32 


C26 




Press 


Lire Condition: 


V =LC 


a am 


< M r’a! 


.> [ rJ ci.ii! 


T emme rat are 


_ /,T) 

- O 


?o 

. . 'j 


< Deg C ) 


N d. t l 


t 1 o w Kate ^ .. ) 


* 30 


. 00 


N<i r . 0 r 


Celoci ty 


= A 


. A 1 


C*i • s ) 


Heat 


F i ax 


~ J • o 


ii.SE 


+05 (W/m 2) 



UNCERTAINTY ANALYSIS: 




VARIABLE 


PERCENT UNCERTAIN 


Mass Flow Rate. Mo 


1 .44 


Reynolds Number. Re 


t 57 


Neat Flux, q 


0.50 


Log- Mean- 'em Di^f. LMTD 


0.43 


Wall Resistance. Rw 


2.67 


Overall H . T . C . , Uo 


.66 


Water-Side H.T.C.. Hi 


1 .49 


Steam-Side H.T.C., Ho 


2.31 
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DATA F GR THE UNCERTAINTY ANALYSIS: 



c j i e name : 




-32 ASS 


Pressure Condition: 




Atmospheric ' ! 0 1 k° 


Steam Temperature 


- 


38.33 C> 


Hater Fiou Rate <«> 


= 


20 . DO 


Water Velocity 


- 


1 . ! A (r,/c.) 


Heat Flux 




6 . 502E + 05 < H/m ’ 2> 


UNCERTAINTY ANALYSIS: 


VARIABLE 




D ERCEN T UNCERTAINTY 


Mass r j ou Rate. K, d 




5.45 


Reynolds Humber, Re 




S . 5C 


Hear Flux, a 




1.74 


i_og-Nean- Ten D 1 F r . i_fi T D 




0.07 


Hall Resistance. Ru 




2.67 


Overall ri . T . C . . Uo 




1 . 74 


Hater “Side H.T.C. , Hi 




4.47 


S' t o an - $ 1 c e H . T . C . . Ho 




or, 77 
• t L_ 


DATA “OR THE UNCERTAINTY 


ANAL YS 


7$' ; 








File Vane : 




u 3 2 h 8 


Pressure Condition: 




ft f nof.oner i c HQt V 


Steam Temperature 


- 


?H ,3? (Deg O 


Hater r i ou K are < 7 > 


= 


8 0 . 0 0 


Hater Velocity 


=■ 


4.^0 ( n / s ) 


Heat Flux 


~ 


1 . 164E+06 ~<U/r 2 2 > 


UNCERTAIN T Y ANALYSIS: 


VARIABLE 




PERCENT UNCERTAINTY 


Mass Fiou Rate, Md 




i .44 


Reynolds Number, Re 




! .59 


Heat Flux, q 




0 . 48 


Log-Mean- Ten Diff. LMTD 




0.14 


Hail Resistance, Ru 




2.67 


Overai 1 H . T . C . . Uo 




50 


Hater -Side H.T.C. . Hi 




1 .50 


Steam-Side H.T.C., Ho 




2.80 



•’a) 
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DA T A -OR THE UNCERTAINTY ANALYSIS: 



Rile Name : 

P r os s ur o Condition: 
Steam Temperature 
Water Flou Rate (*/<) 
Water Velocity 
Heat Fiu.x 



UNCERTAINTY ANALYSIS: 
VARIABLE 

Mass Fiou Rate. Md 
Reynold?. Number, Re 
heat Flux, c; 

Log-Mean- Ten Diff. L.MTD 
Wall Resistance. Ru 
Overall H . T . C . . U o 
Water -Side H.T.C.. Hi 
Steam-Side H.T.C.. Ho 



F39V83 

Vacuum < 1 1 k Pa > 

= 48 50 (Deo O 

= 20.00 

= l 03 ( m / <= ) 

= 2 . 579E+05 (W/m 2) 



D ERCEN1 UNCERTAINTY 

5.42 
5 . -5 
1 . 73 
0.23 
A. 39 
I .75 

4.43 
13.85 



DATA FOR THE UNCERTAINTY 


ANALYSIS: 




file : iane : 


F3SV83 




"Tessure Condition: 


Vacuum 


( ! ! ORa) 


Steam Temperature 


= 48.4! 


< Deg C> 


Water Flou Rate <7) 


= 80.00 




Hater Veiocir. v 


_ *1 r L 

I' . ?>L 


( n / t > 


neat F l ux 


= 4.440E 


0 5 ( W / m ’ 2 > 



UNCERTAINTY ANALYSIS : 




VARIABLE 


PERCENT UNCERTAINTY 


Has? Flou Rate, Hd 


: .43 


Reynolds Number, Re 


! . j5 


Heat Flux, q 


0.50 


Log-Mean- Tem Diff. LI1TD 


0 .48 


Wall Resistance. Ru 


4.39 


Overal 1 H . T . C . , ! Jo 


.69 


Water -Side H.T.C.. Hi 


1 .47 


Steam-Side H.T.C., Ho 


1 .89 
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DATA FOR THE UNCERTAINTY 


ANALYSIS: 


File Name : 

Pressure Condition: 
Steam Temperature 
Water Flow Rate <*/> 
Water Velocity 
Heat Flux 


F39A98 

Atmospheric <!C kP?.) 

= A 9 . 9 I ( [) P g C ) 

= 20 . 00 
= ! . 0 A ( m / s > 

- 8. 13 IE + 05 < W / m ' 2) 


UNCERTAINTY ANALYSIS: 
VARIABLE 


D ERCEN T !JNCEP T AINTY 


Mass Flow Rate. Md 
Reynolds Number. Re 
Heat Flux, q 
Log-Mean- '-a Dif f . LMTD 
Wall Resistance. Rw 
Ove ra i 1 H . T . C . . Uo 


5 . 29 
5. A3 
1 . 72 
J.07 
A . 39 
1 .72 


Water-Side H.T.C.. Hi 
St earn- Side W.T.C.. Ho 


A . 42 

:s.s? 


DATA FOR THE UNCERTAINTY 


ANALYSTS: 


File Name : 

Pressure Condition: 
Steam Temperature 
Water ~ 1 oui Rate (’/,'* 
Water Velocity 
Heat Flux 


P39A98 

Atmcspneric <!0I kPa) 
-- 99.78 ( Deo C) 

= 80.00 
= 3 , j4 ( m/e ) 

= . 4 0 9 E + 0 n ( W / m 2 ) 


UNCERTAIN' i ANALYSIS: 
VARIABLE 


■PERCENT UNCERTAINTY 


Mas s Flow mate, M d 
Reynolds Numoer, Re 
Heat Flux, q 
Log-Nean-Tem Diff. LHTD 
Wall Resistance, Rw 
Clve ral 1 H . T . C . , Uo 
Water-Side H.T.C., Hi 
Steam-Side H.T.C.. Ho 


1 .42 
1 .54 
0.48 
0. 15 
4.39 
.50 
1 . 4b 
1 .89 



141 



LIST OF REFERENCES 



Yau., K. K . , Cooper, J. R. , and Rose, J. W. , "Effects 
of Fin Spacing and Drainage Strips on the Condensation 
Heat-Transfer Performance of Horizontal Low 
Integral-Fin Tubes". Fundamentals of Phase C hange ; 
B oil ing and Condensation, HTU-VoI 38, C. T. AvedTsian 
an3“T. M. Hudy (Eds.f EBME, 1984. 



Wanniarachchi, A. S. , Marto, P. J., and Rose, J. W. , 
"Filmwise Condensation of Steam on Externally-Finned 
Horizontal Tubes", Fundamentals of Phase Change: 
B oil ing and Condensation, HTF = 75I. 3B, C. T. Avedisian 
Ind"T7 M7“Pudy~lEds7r _ IS?iE, 19 84. 



Honda, H., Nozu, S., Mitsumori, K. , "Augmentation of 
Condensation on Horizontal Finned Tubes By Attaching 
Porous Drainage Plates", Proc. AS ME- JS ME Th erma l 

E ngine eri ng Conference, Hawaii, 1983. 



Krohn, R. L., An Exper imental A ppara tus to Stud y 
Enhanced Conden sat ion He at- Tr ansf er of Steam on 

Hor izon tal TuFes, IT. S. Tfesis, faval^Postgraduate 
Befool, Monterey, California, June, i982. 



Graber, K. A. , Condensation Heat Transfer of Steam on 
a Single Horizontal Tufe, M. B. Tfesis, Naval 
Postgraduate Befool, Monterey, California, June, 1983. 



Poole, F7 . M., Filmw ise Condensation of Ste am on 

Exter na lly- Fin ned Horizontal T ufes . ff7B. Tfesis, TIaval 
Postgraduat “Befool, Honterey, California, December, 
1 9 83 • 



Georgiadis, I. V., Filmwise Condensation of Steam on 
Low Integr al - fin ned Tufes, HTB. Tfesis, Paval 
Postgraduate BcfooTJ Monterey, California, September, 



Katz, D. L. , Hope, R. E. , and Dasko, S. C. , Liquid 
Rete n ti on on F inned Tube s, Dept. of Eng. Researcf, 
Oniverslty of Michigan , fnn Arbor, Michigan, Project M 
592, 1946: 



Gregorig. R. , "Hau tkondensation an Feingewellten 
Oberflachen bei Berucksichtigung der 
Oberfachenspannungen" , Zeitschrift fur Angewantde 
Mathematik und Physic , Tol. 7, T954, pp. 36-49. 
Translation Fy D. K. Edwards. 



F.udy, T. M., and Webb, R. L., "Condensate Retention of 
Horizontal Integral-Fin Tubes, Advances in Enhanced 
Heat-Transfer", 1981, HTD-Vol. Iff, Presented at 2Ufn 
2UTE FTaffional Heat-Transfer Conference, Milwaukee, 
Wisconsin, August, 1981. 



Rudy, T. M., and Webb, R. L. , "An Analytical Model to 
Predict Condensate Retention on Horizontal. 
Integral-Fin Tubes", ASME/JSME Thermal Engng. Joint 
Conf . , Vol. 1, Marcn 20^2ET, T9ff3. 



Rifert, V. G. , "Steam Condensation on Profiled 
Surfaces". Heat and Ma ss- Transfer Processes in Porou s 
Media With Phase Transformation, Academy of Science, 

ffssrr ir - ff.~ryicov - TEaT 7 - Miisic 7 ~i 9 82, pp. 1 49-170. 



Owen, R. G., Sardesai, R. G., Smith, R. A., and Lee W. 
C . , Gr avity Co ntr olled Condensation on a Low-Fin Tube, 
I.Chem.E. Symposium Series Ho. 75. 



Rudy, T. M., and Webb, R. L. , "Condensate Retention on 
Horizontal Integral-Fm Tubes", ASME Journal of Heat 
Trans fer. Vol. 107, 1985. 



Beatty, B. 0., and Katz, D. L., "Condensation of 
Vapors on Outside of Finned Tubes", Ch emica l 
E ngineer ing Pr ogress . Vol. 44, No. 1, January, ~ TP 48 . 



Karkhu, V. A., and Borovkov, V. P. , "Film Condensation 
of Vapor at Finely-Finned Horizontal Tubes", Heat 
Transfer-Soviet Research, Vol. 3, No. 2, March-Apnl 
T97T7 



Zozulya, N. V., Karkhu, V. A., and 
"An Analytic and Experimental Study 
in Condensation of Vapor on Finned 
Trans fer -Soviet Resea rch . Vol. 9, No 



Borovkov , V . P. , 
of Heat Transfer 
Surfaces", Heat 
. 2, March-Aprll 



Nader, W. K. , "Extended Surface Heat Transfer With 
Condensation", He at T ransfer 197 8 , Vol. 2, August, 



Patankar, S. V., and Sparrow, E. M., "Condensation on 
an Extended Surface", Transactions of the ASME, Vol. 
101, August, 1979. 



Rifert, V. G. , "A New Method for Calculating Rates of 
Condensation on Finned Tubes", Heat Tra nsfer-Sov iet 
Resea rch . Vol. 12, No. 3, May- June, 198(17 



Rudy, T. M. , and Webb, R. L. , "Theoretical Model for 
Condensation on Horizontal, Integral-Fin Tubes", Heat 
Transfer, Seattle, AIChE Symp. Ser., Vol. 79, No. 225, 



Webb, P. . L. , Keswani, S. T. , Rudy, T. M. , 
"Investigation of Surface-Tension and Gravity Effects 
in Film Condensation", P roceeding s of 7th 
I nt er nat io nal Heat-Tr ans fer Conference, Hunich, Fed. 
Hep. of Germany, Sep"f . d-TOT TS82, Hemisphere 
Publishing Co., Washington D.C., Vol 5, pp. 175-1 80. 



Rudy, T. M. , A Th eo r etica l and Experimental Study of 
Co nden s at ion on Single, Integral- Tin luges, PK.I5. 
Thesis. Eeparfmenf of Mechanical "Engineer mg. The 
Pennsylvania State University, UNiversity Park, Pa. , 
May, 1982. 



Honda, H. , and Nozu, S., A Prediction M et hod for 
Heat-Transfer Purring Film Condensation in Horiz onta l 
low T nte gr a I-Fm~ TuE es , Submitted" for presentation at 
the A" SHE Winter Annual Meeting, New Orleans, December, 
1984. 



Webb, R. L. , Rudy, T. M. , and Kedzierski, M. A., 
"Prediction of the Condensation Coefficient on 
Horizontal Integral-Fin Tubes", Journal of Heat 
Transfer, Vol. 105, 1985. 



Rudy, T. M. , Kedzierski, M. A., and Webb, R. L. , 
"Investigation of Integral-f in-Ty pe Condenser Tubes 
for Process Industry Applications", First U. K. 
Nati o na l Conference on He at Transfer, The Institution 
of Chemical Engineers Symposium Series No. 86, 1984. 



Adamek, T., "Bestimmung der Kondensationgrossen auf 
feingewellten Oberflachen zur Auslegung aptimaler 
WandjDrof ile" ^ Warme-u n d-Sto f f uh ert raqunq , vol. 15, pp 



Incorpera, F. P. , and DeWitt, D. P., Fundamentals of 
H eat Transfer, Wiley, pp 406-407, 1981. 



Fujii, T. , Honda, H., "Condensation of Steam on a 
Horizontal Tube", Condensation Heat Transfer, ASMS, 
New York, 1979. 



Webb, R. L., Haman, L. L. , and Hui, T. S., "Enhanced 
Tubes in Electric Utility Steam Condensers", Heat 
Transfer in Heat Rejection Systems, HTD-Vol. 37, ~S. 
Sen^upFa and Y7 HussalTi lEds.f/ ASMS, 1984, pp. 



Fujii, T., Wang, W. , Koyama, S. , and Shimizu, Y. , 
H ea t-T ra ns f er Enha nce ment for Grav it y Contr aile d 
Co ndensa f Ion on a Hori zon t al Tube hy Coiled W ires , 
Beigmg Conference, TSE5. 



Kline. S. J. , and McClintock, F. A., "Describing 
Uncertainties in Single-Sample Experiments", Mech . 
Eng. , Vol. 74, pp 3-8, January 1953. 



BIBLIOGRAPHY 



Katz, D. L. , and Geist, J. M . , "Condensation on Six Finned 
Cubes in a Vertical Row", Transactio ns of the AS MS, 
November, 1948. 

Mori, Y. , Higikata, K. , Hirasawa, S. . and Nakayama, W. , 
"Optimized Performance of Condensers With Outside Condensing 
Surfaces", Transactions of the ASHE, Vol. 103, February, 
1981. 

Pearson, J. F. , and Withers, J. G. , "New Finned Tube 
Configuration Improves Refrigerant Condensing", ASHRAE 
Journal , June, 1969. 

Shklover, G. G., Mil’Man, 0. 0., Baskov, V. S. t Ankudinov, 

G. A., "Heat Transfer in Condensation or Steam on 
Finely-Finned Horizontal Tubes", Heat Transfer-Soviet 
Re se arch , Vol. 13, No. 2, March-April T^ffT. 



145 



INITIAL DISTRIBUTION LIST 



No. Copies 

1. Defense Technical Information Center 2 

Cameron Station 

Alexandria, Virginia 22314 

2. Library, Code 0 142 2 

Naval Postgraduate School 

Monterey, California 93943 

3. Department Chairnan, Code 69 1 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey, California 939 43 

4. Professor P. J. Marto, Code 69Mx 5 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey, California 93943 

5. Dr. John W. Rose 1 

Department of Mechanical Engineering 

8 ueen Mary College 
niversity of London 
London El 4NS 
England 

6. Dr. A. S. Fanniarachchi, Code 69Wa 1 

Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, California 93943 

7. Dr. Win Aung 2 

Program Director for Heat Transfer 
Division of Engineering 
National Science Foundation 
Washington, D.C. 20008 

8. Lt. Frederick A. Flook, USN 2 

Long Beach Naval Shipyard 
Long Beach, California 90822 

9. It. Ioannis Georqiadis, H. N. 1 

67. Peloponisou St, Agia Paraskevi 
Attiki, Athens 
Greece 

10. Lt. Evangelos Mitrou, H. N. 1 

SMC #2777 

Naval Postgraduate School 
Monterey, California 93943 

11. Mr. Krohn, R. 1 

Supervisor of Configuration Management 
Clinton Power Station 
Decatur, Illinois 62526 



146 



1 

s 



v 

) 



\ 

* 






1 3 5 3 7 



5 











I 














:i2t 



o 



Thesis 

„F522U3 Flook 

c.l Film condensation 

of steam on externally 
finned horizontal 
tubes. 



/ 



?: 



Thesis 

F52243 Flook 

c.l Film condensation 

of steam on externally 
finned horizontal 
tubes. 



